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CHAPTER  1 

IMTBQPUCTIQH 


This  report  is,  in  essence,  concerned  with  scheduling 
theory.  The  concern  takes  two  forms.  First,  there  is  a 
practical  engineering  scheduling  problem  that  needs  to  be 
solved.  To  study  this  problem  we  have  drawn  upon  a new 
theoretical  approach  to  scheduling.  And,  this  is  our  second 
concern,  the  development  of  a theory  that  indeed  has 
practical  implications. 

The  problem  we  address  is  that  of  scheduling 
dlal-a-ride  transportation  systems.  A solution  would  be  the 
development  of  a methodology  for  analysing  and  designing  the 
scheduling  algorithms.  Whilst  aiming  at  a solution,  we 
investigate  to  what  extent  a recent  asymptotic  probabilistic 
technique  for  the  solving  of  hard  combinatorial  optimisation 
problems  is  of  real  Interest.  This  requires  that  a 
matbematloal  result  be  generalized  for  a number  of  problem 
formulations. 
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1 . 1 The  Dlal-a-Hide  Problem 


During  the  past  decade  there  has  been  some  interest  in 
the  planning  of  innovative  public  transportation  systems. 
An  area  which  has  received  substantial  attention  is  that  of 
'demand  responsive  transportation*  [19]. 

One  of  the  outcomes  of  this  research  has  been  the 
' dial-a-ride ' proposal.  A dlal-a-ride  transportation  system 
is  somewhere  in  the  range  between  a rigid  bus  system  and  a 
flexible  taxicab  system,  and  ideally  provides  large  numbers 
of  passengers  with  personalized  service.  Passengers  request 


service  - to 

be  taken 

from  an 

origin  to 

a destl 

nation 

- by 

telephone.  At 

this 

time 

of  req 

uest , 

an  e 

stimate 

of  the 

time 

of  collection 

and 

time 

of  del 

iver  y 

is  q 

uoted . 

Small 

buses 

travel  about  the  region,  collecting  and  delivering 
passengers;  during  his  trip,  a particular  passenger  may  be 
transferred  to  another  bus.  At  any  time  a bus  may  have  many 
passengers  aboard. 


The  central  mathematical  problem 
transportation  system  is  one  of  schedulin 
required  that  will  decide  to  which  bus 
buses)  a particular  passenger  should  b 
his  trip  should  take  place.  A number  o 
been  proposed  to  do  this  scheduling, 
computer  aided;  the  computer  is  used 
on-line  decision  making  or  as 
evaluating  the  system  'state'. 


in  a dial-a-ride 
g:  an  algorithm  is 

(or  sequence  of 
e assigned  and  when 
f algorithms  have 
Host  of  these  are 
for  the  detailed 


an  information 


source 
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We  shall  be  concerned  with  dial-a-rlde  scheduling 
algorithms  at  a fundamental  level,  avoiding  many  of  the 
modifications  and  adjustments  that  must,  of  course,  be  made 
in  an  actual  implementation  (we  shall  discuss  some  of  the 
practicalities  in  Chapter  IV).  There  are  two  fundamentally 
very  different  approaches  to  scheduling  which  we  shall 
consider.  Let  us  describe  the  skeletons  of  these 
approaches . 

The  first  is  a system  developed  at  M.I.T.  by  Wilson 
et.  *1.  [21,22];  experimental  versions  have  been 

ented  since  1972,  and  currently  one  is  being  tested  at 
ester.  New  York.  The  underlying  trait  of  this  scheme  is 
its  search  procedure  for  allocating  passengers  to  buses.  At 
each  instant  of  time  each  bus  has  associated  with  it  a 
'prospective  route',  given  by  an  ordered  sequence  of  future 
stops  (which  may  be  origins  or  destinations ^ ^ ^ , and  an 
estimated  time-of-arrlval  at  each  of  these  stops.  With  each 
stop  too  there  is  a 'latest  time-of-arrival ' - the  time 
which  has  been  quoted  to  the  associated  passenger.  A new 
Incoming  passenger  must  be  allocated  to  a bus.  This 
Involves  inserting  his  origin  and  destination  into  the 
prospective  route  of  one  of  the  buses.  For  each  new  demand, 
the  scheduler  searches  through  every  possible  insertion  on 
every  one  of  the  buses  and  chooses  the  best,  so  as  to 

(1)  This  scheme  requires  no  transfers  for  passengers:  each  t 

passenger  must  be  delivered  by  the  same  bus  that  collects  ■ 

him.  I 
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( 2 ) 

minimize  a certain  criterion  function  • Basic  to  the 
design  of  algorithms  here  is  the  choice  of  the  criterion 
function,  the  method  of  quotation,  and  the  travel-time 
prediction.  It  appears  that  the  procedure  is  very  sensitive 
to  these  choices,  as  well  as  to  variations  in  the  parameters 
[21,22,18]. 


At  the  other  extreme  is  a very  different  approach 
towards  the  scheduling.  It  is  described  by  the  following, 
which  we  shall  refer  to  as  the  'Michigan  Scheme'  since  a 
related  system  has  been  evolving  at  Ann  Arbor,  Michigan 
[17].  The  region  in  which  the  system  operates  is 
partitioned  into  a number  of  subregions,  r ^ , r2 , . . . rn . In 
each  region  rj  there  is  a bus  that  travels  only  in  r^.  A 


{2)  The  criterion  is  a function  that  concerns  both  present 
and  future  passengers.  Present  passengers  are  interested  in 
their  wait  time,  w,  ride-time,  r,  total  travel  time,  w-«-r, 
pickup-time  deviation.  Dp,  and  delivery-time  deviation,  Dd 
(these  are  deviations  from  the  quoted  times).  Future 
passengers  are  acknowledged  by  taking  into  consideration  the 
increase  in  tour-length,  DT.  Thus,  for  each  possible 
assignment  one  can  evaluate 

a^flCw)  ♦ a2f2(<^)  ♦ a3f3(w+r)  ♦ a4fi((Dp)  a5f5(Dd)  (1) 

for  every  current  passenger,  and  a similar  value  for  the 
passenger  being  assigned.  Here,  the  ai'*’  weighting 

parameters  and  the  ^i ' * s**®  functions  (generally, 

fj[(z)  s z or  x^  for  all  1).  Summing  (1)  for  all  the 
passengers  and  adding  also 
bDT 

where  b is  another  weight,  yields  the  criterion  function  to 
be  minimized. 

The  time  to  a point  P which  is  quoted  a passenger  is 
determined  by 

c.E(tp)  <4-  d 

where  E(tp)  is  the  expected  direct  travel-time  to  P and  c 
and  d are  fixed  parameters.  In  some  versions  of  the 
algorithm  the  quoted  times  are  retained  as  hard  constraints, 
which  must  never  be  violated  by  subsequent  assignments. 
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larger  'line-haul'  bus  connects  the  regions  by  travelling 
along  a fixed  route,  stopping  at  certain  'transfer  points', 
PliP2****Pm  each  region.  This  is  shown  schematically 

in  Figure  1.1. 


Figure  1.1  The  "Michigan  Scheme" 

A passenger  requiring  to  go  from  region  r^  to  region  rj 
is  collected  by  the  bus  in  rj^,  transfers  onto  the 
line-haul  bus  at  p^^,  alights  at  pj  and  is  delivered  by 
the  bus  in  rj.  The  regional  buses  visit  their  transfer 

I 

• points  every  15  minutes  (say).  The  regions  are  small  enough 

t 

I that  drivers  can  quite  easily  choose  an  optimal  or 

near-optimal  path  between  visits  to  their  transfer  point. 
The  times  which  are  quoted  to  passengers  for  collection  and 
delivery  are  simple  estimates  of  travel  time  given  the 
load^^\  Under  a heavily  loaded  system  a passenger  may  have 
to  wait  his  turn  before  being  collected,  as  the  regional 


(3)  These  quotations  are  not  used  in  the  scheduling 
procedure . 
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buses  can  make  only  a limited  number  of  stops  between  visits 
to  the  transfer  point. 

We  do  not  intend  to  give  here  a detailed  evaluation  of 
present  dial-a-ride  algorithms  - the  reader  is  referred  to 
[10]  - but  some  critical  comments  are  in  order.  These  focus 
on  three  aspects:  the  cost,  the  optimality  and  the 
practical  performance.  We  base  our  initial  observations  on 
the  Rochester  approach,  since  there  has  been  a great  deal  of 
research  into  this  and  related  schemes. 

Basic  to  the  Rochester  approach  is  a long  and  hard 
search.  Typically,  the  computational  effort  required  for 
such  searches  grows  exponentially  with  the  problem  size  (see 
the  comments  in  Section  2.1),  and  we  should  expect  it  to  be 
very  expensive. 

A scheme  might  be  Justified  if  a certain  performance 
level  can  be  guaranteed.  The  Rochester  search  is  a local 
one  [5],  both  with  respect  to  time  and  with  respect  to  the 
passengers.  That  is,  only  a single  passenger  is  assigned 
(or  at  best  only  a small  set  of  passengers  is  reassigned) 
when  the  schedule  is  updated.  Further,  optimality  (with 
respect  to  the  criterion  function)  is  ensured  for  only  that 
point  in  time.  For  local  searches  it  is  very  rare  that  good 
performance  can  be  guaranteed.  Recent  research  in  [22]  has 
focused  on  the  criterion  function.  However,  there  is  no 
guarantee  that  the  performance  of  the  system  as  a whole 
(measured,  say,  in  terms  of  average  travel-time  and/or 
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( j| ) 

variance  of  travel-time  ) will  be  improved  If  the  utility 
functions  of  all  current  passengers  are  maximized  whenever 
decisions  are  made.  This  Is  true  even  If  reasonably 
acceptable  utility  functions  for  the  passengers  could  be 
represented,  Itself  a notoriously  messy  problem. 

Lastly,  what  about  Its  practical  performance?  In  the 
early  stages  of  Implementation,  particularly  when  the  system 
was  heavily  loaded,  there  were  often  roundabout  devious 
routes  for  passengers,  with  resulting  customer 
dissatisfaction  [16].  In  general,  it  appears  that  the 
results  on  performance  are  Inconclusive. 

In  contrast,  the  simplicity  of  the  scheme  at  Michigan 
is  very  attractive.  The  scheduling  cost  is  relatively  low, 
since  the  computer  Is  used  only  to  store  the  demands  and 
access  them  efficiently.  (Optimal  tours  for  the  regional 
buses  might  also  be  determined  on-line  by  a computer.  Then 
each  tour  Is  but  a small,  simple  'travelling  salesman  tour* 
- see  Chapter  II  - where  the  distances  between  points 
satisfy  the  triangle  Inequality.  For  these  problems  there 
are  some  recent  efficient  heuristic  algorithms  [12].) 

In  practice  the  Michigan  scheme  has  been  found  to  work 
well  [17].  But  there  are  still  some  vague  questions.  Can 
one  Justify  the  use  of  the  Michigan  approach  with  respect  to 
some  optimality  criterion;  and  how  should  one  design  the 
regions,  the  fixed  routes  and  the  timing  of  transfers? 


(4)  These  were  the  measures  loosely  used  In  [22]. 
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SiilX.  Approach 

Our  investigation  will  be  an  analytic  one  into  some  of 
the  basic  aspects  of  dial-a-ride  algorichms.  We  shall 
derive  a class  of  algorithms  - or,  more  correctly,  an 
approach  towards  the  design  of  algorithms  - for  which 
optimality  can  be  measured  in  a certain  precise  sense.  This 
sense  is  asymptotic  (in  the  number  of  passengers  who  require 
service)  and  probabilistic  (so  that  the  probability  is  high 
that  on  any  given  day  the  algorithm  will  perform  well). 

The  approach  taken  is  one  of  preplanning  at  a global 
level.  Very  loosely,  it  is  based  on  the  following 
principle.  Even  though  each  passenger  is  unique,  with  his 
own  required  origin,  destination  and  time-of-dellvery , in  a 
large  system,  where  there  are  a large  number  of  passengers, 
we  can  predict  quite  closely  the  behaviour  patterns  of  the 
set  of  all  passengers.  This  is  the  'equalizing'  effect  of 
the  law  of  large  numbers  that  has  been  observed  in  many 
physical  phenomena,  the  classic  example  being  the 
thermodynamic  principles  for  the  behaviour  of  gases. 

To  analyse  the  problem  theoretically  it  is  necessary  to 
abstract  the  essentials  and  to  consider  an  idealized 
version.  Chapters  II  and  III  present  an  analysis  of  models 
of  the  dial-a>rlde  problem  that  focus  upon  its  combinatorial 
nature.  They  study  separately  static  versions  - in  which 
demands  for  service  are  all  available  at  the  start  - and 


dynamic  versions  - in  wbiob  the  demands  arise  over  time 


I 

t 


¥ 
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For  the  nodela,  we  obtain  'asyaptotlcally  optiaal' 
algoritbaa  that  minimize  simple  distance  or  average 
flow-time  criteria,  and  we  evaluate  suboptlmal  schemes. 

Despite  the  idealization,  many  qualitative  insights 
result,  and  in  Chapter  IV,  we  return  to  the  real  problem. 

Here,  drawing  upon  the  theoretical  results,  we  obtain  their 
implications  for  the  real  problem  and  propose  an  approach 
towards  the  design  of  dial-a-ride  systems.  The  approach  is 
of  interest  then,  because  it  can  be  theoretically  Justified, 
and  also  because  it  has  many  attractive  practical  features. 
Furthermore,  the  techniques  developed  in  Chapters  II  and  III 
provide  us  with  a powerful  analytical  tool  for  use  in  the 
design  process.  We  can  Investigate  changes  in  performance 
when  parameters  (e.g.  the  number  of  buses,  the  size  of  the 
region  and  many  others)  are  varied.  Proposed  schemes  can  be 
easily  evaluated  without  the  need  to  resort,  at  this  basic 
level,  to  simulation. 

It  is  interesting  to  note,  at  this  point,  that  the 
'Michigan  Scheme'  described  above  can  be  considered  as 
belonging  to  our  class  of  algorithms.  Thus,  we  are  able  to 
say  in  what  sense  it  is  optimal  as  well  as  to  gain  some 
valuable  insights  towards  its  improvement.  The  attractive 
simplicity  of  the  Michigan  scheme  is  common  to  all  of  our 
proposed  algorithms.  ^ 

I 


CHAPTER  II 


STATIC  MS  PROBLEMS 


In  this  chapter  we  begin  a theoretical  analysis  of 
dial-a-ride  bus  systems.  In  Section  2.0,  we  introduce  the 
idealized  'bus  problems'  which  we  shall  be  developing,  and 
state  the  basic  assumptions  and  terminology.  Our  approach 
is  based  on  an  interesting  theorem  of  Beardwood,  Halton  and 
Hammersley  [2]  and  on  a recent  paper  of  Karp  [9];  we 
describe  these  results  within  the  context  of  combinatorial 
optimization  problems  in  Section  2.1.  Sections  2.2  and  2.3 
present  a theory  for  single-  and  multiple-bus  problems. 

Throughout  this  chapter  the  presentation  is  informal, 
with  only  heuristic  Justifications  and  proofs.  Detailed 
proofs  of  the  results  are  delegated  to  an  appendix. 


IdealiMd  Ma  Problems 

The  travelling  salesman  problem  (TSP)  - see  the 
definition  below  - has  received  much  attention  in  the 
operatlons-research  literature.  Not  least  among  the  reasons 
for  this  is  the  fact  that  the  TSP  can  be  regarded  as  a 
prototype  for  many  realistic  problems.  Clearly,  there  is  a 
nontrivial  relationship  between  the  TSP  and  the  dial-a-rlde 
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problem,  and  in  tbla  chapter  we  shall  exploit  this 
relationship.  First  let  us  describe  the  problems  to  be 
considered,  and  state  the  basic  assumptions  and  our 
terminology. 


We  are  given  a bounded  region  R in  the  plane,  with 
area  a.  Demands  arise  in  R.  A demand  p s (o,d)  is  a 
pair  of  points  in  R with  origin  o and  destination  d.  A 
set  of  m buses  Is  at  our  disposal  for  meeting  these 
demands  - l.e.  for  visiting  the  set  of  points.  A tour  for 
a bus  is  the  sequence  of  points  which  it  visits.  Before  a 
bus  can  visit  a destination  the  relevant  passenger  must  be 
on  board,  and  bo  the  corresponding  origin  must  already  have 
been  visited.  We  refer  to  this  as  the  'feasibility 
constraint';  tours  which  satisfy  it  are  termed  'feasible'. 


We  distinguish  between  static  and  dynamic  versions  of 
the  problem.  In  the  static  version  we  are  given,  at  time 


tsO , a 

collection  of  n 

demands 

and 

we 

are 

to  devise 

feasible 

tours  for  the 

buses  so 

that 

all 

2n 

points  are 

visited.  In  dynamic  problems,  to  be  discussed  in  Chapter 
III,  the  demands  arise  as  time  progresses  according  to  some 
random  process.  Again  feasible  tours  are  to  be  devised  (of 
course  a point  can  be  visited  only  after  the  demand  arises). 
The  problem  faced  is  to  devise  tours  .10  that  a certain 
criterion  is  minimized.  The  criterion  will  usually  be  a 
function  of  tine  or  of  distance  travelled. 


r 
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We  assuoe  that  the  buaes  are  all  of  infinite  size.  The 
distance  between  two  points  in  R is  the  euclidean  length 
of  the  straight  line  Joining  then,  so  if  R is  not  convex 
the  tour  might  leave  R.  The  buses  travel  at  unit  speed, 
and  no  time  is  wasted  when  a passenger  embarks  or  alights. 
Also,  transfers  are  instantaneous  (see  Section  2.3).  These 
assumptions  will  be  discussed  further  in  Chapter  IV. 

An  instance  of  a static  bus  problem  of  size  n is 
specified  by  a set  of  n demand  pairs.  Our  approach  is 
probabilistic  and  we  must  define  a probability  distribution 
from  which  the  problem  Instances  are  drawn.  We  assume  (for 
simplicity,  although  this  can  be  considerably  weakened)  that 
all  origins  and  destinations  are  drawn  Independently  from 
the  uniform  probability  distribution  over  R.  We  refer 
loosely  to  this  fact  by  saying  that  a particular  problem 
instance  (of  size  n)  is  'random'. 


We  shall  also  be  making  use  of  the  well-known 
travelling  salesman  problem.  An  instance  of  this  problem 
(the  euclidean  version)  is  given  by  a set  of  n points 
within  R.  We  are  required  to  find  a patb^^^  which  passes 
through  these  n points  and  which  has  shortest  length.  In 
a 'random'  instance  of  the  TSP,  the  n points  are  drawn 
independently  from  a uniform  probability  distribution  over 
I. 
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► 

In  the  algorithms  to  be  considered,  the  region  R is 
partitioned  into  subregions  (this  emphasis  will  be 


Justified).  When  we  say  a bus  'visits'  certain  points  in  a 
subregion  we  mean  that  it  enters  the  subregion  and  performs 
an  optimal  travelllr  salesman  tour  on  the  designated  points 
there.  When  the  ^ u.*.  >,hen  visits  points  in  another  region, 
it  travels  to  the  closest  of  the  new  points  and  similarly 
performs  an  optimal  travelling  salesman  tour  in  that  region. 


If  we  think  of  R as  being  divided  into  m equal 
subregions  each  of  area  a/m,  then  any  given  demand  pair 
originates  in  any  particular  region  with  probability  1/m, 
and  has  destination  in  any  region  with  probability  l/m.  If 
there  are  n demands,  with  n a large  number,  then  with 
high  probability  there  are  approximately  n/m  origins  and 
n/m  destinations  in  each  region.  This  follows  from  the 
strong  law  of  large  numbers  [14].  However,  we  might  loosely 


assume  that  there  are  exactly  n/m  origins  and  n/m 

destinations  in  each  region.  This  is  an  example  of  the  sort 
of  imprecision  that  exists  in  the  following  pages.  Rigorous 
statements  and  detailed  proofs  of  the  results  are  collected 
in  the  appendix,  so  we  are  free  to  concentrate  here  on  the 
spirit  of  the  ideas.  It  should  be  kept  in  mind  throughout 
that  inexactitudes  as  that  above  hold  only  asymptotioally 
(in  n)  with  probability  1. 
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Our  analysis  will  begin  with  a very  simple  bus  problem, 
and  the  results  will  then  be  extended  to  more  complex 
problems  with  corresponding  Increased  realism.  First, 
however,  we  digress  to  describe  the  aspects  of  scheduling 
theory  upon  which  our  approach  is  based. 


2. 1 Combinatorial  OptlBization  AJil  Ikfi  IraYftZllng  Salaaaaii 
Probiaa 


Recently  there  have  been  some  theoretical  advances  in 
the  understanding  of  combinatorial  optimization  problems. 
Most  Important  has  been  the  acknowledgement  that  there 
exists  a class  of  'hard  problems'  (NP-complete  problems). 
These  are  hard  in  the  sense  that  the  computing  time  required 
to  find  their  solution  by  any  known  algorithm  explodes 
exponentially  as  the  size  of  the  problem  increases.  For 
these  hard  problems,  then,  there  seems  to  be  no  way  to 
avoid,  essentially,  the  enumeration  of  a very  large  number 
of  possible  alternatives.  This  is  inefficient;  an 
'efficient'  algorithm  would  be  one  with  an  execution  time 
behaving  as  p(n),  a polynomial  function  of  n.  These 
notions  have  been  made  mathematically  precise  • see  iho 
et.  al.  [1],  Coffman  [5]. 

For  example,  the  TSP  has  been  shown  to  be  NP-oomplete 
116],  and  there  are  no  known  efficient  algorithms  for  its 
solution.  To  solve  it,  we  need  to  evaluate  very  many  of  the 
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n/2 I possible  orderings  of  the  points,  or  tours. 
Furthermore,  the  existence  of  an  efficient  algorithm  for  the 
TSP  would  imply,  and  be  Implied  by,  the  existence  of 
efficient  algorithms  for  a whole  class  of  hard  scheduling 
and  other  combinatorial  optimization  problems.  At  present, 
this  appears  extremely  unlikely. 

This  Is  indeed  bad  news,  but  practical  problems  must 
still  be  solved.  To  this  end,  note  that  we  required  above 
that  the  algorithm  guarantee  the  optimal  solution.  An 
easier  question  which  might  be  asked  would  be  one  with  the 
optimality  requirement  ‘relaxed.  So,  find  an  efficient 
algorithm  A (if  one  exlstsl)  that  will  provide  a possibly 
nonoptimal  cost  - call  It  c(A)  - but  one  that  Is  close  to 
the  optimal  cost,  c*.  More  precisely,  an  algorithm  A Is 
said  to  solve  the  problem  to  within  the  ratio  y ( y Is  a 
real  number,  larger  than  1}  if 

c(A)  i yc»  (2.1) 

(We  shall  also  use  the  terminology:  A is  c-optlmal  If 

c(A)  1 (1  + e)c*  ) 

This  approach  has  been  found  to  be  very  useful  for 
certain  special  problems,  and  of  no  use  at  all  for  certain 
others.  For  the  travelling  salesman  problem  the  best  that 
has  been  achieved  in  this  sense  is  an  algorithm 
(Christofldes  [3])  that  solves  the  problem  in  polynomial 
time  to  within  3/2. 
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Thls  ratio  is  still  a little  high  and  one  mlgnt  expect 
to  do  better.  For,  given  an  Instance  of  the  travelling 
salesman  problem,  it  is  not  too  bard  to  draw  a path  which 
looks  pretty  much  the  shortest,  and  one  might  then  be 
tempted  to  say  that  with  high  probability  it  is  close  to  the 

optimum.  This  implies  that  we  should  relax  not  only  the 

1 

optimality  requirement  but  also  the  guarantee  which  the 
algorithm  promises.  This  notion  coincides  with  a very 
recent  approach,  suggested  by  Karp  [9]<  Suppose  that  the 
problem  Instance  is  derived  from  a certain  known  probability 
distribution.  Can  we  then  produce  an  efficient  algorithm 
that  will  perform  well  (l.e.  to  within  an  acceptable  ratio) 
with  high  probability?  This  probabilistic  approach  can  best 
be  Illustrated  via  the  TSF. 

Let  us  suppose  that  a particular  problem  instance  of 
the  TSP  is  chosen  by  drawing  n points  independently  from  a 
uniform  probability  distribution  over  the  planar  region  R. 
Let  Lq  be  the  length  of  the  shortest  path  through  these 
n points.  Then  Ln  is  a random  variable. 

( 2 ) 

The  following  theorem  is  due  to  Beardwood 

et.  al.  [2]. 


(2)  The  problem  studied  in  [2]  is  actually  the  eloaed  TSP  in 
wbiob  s closed  tour  is  sought.  It  is  easy  to  see  that  the 
theorem  still  bolds  for  the  open  TSP. 
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Ih^orea  1 

If  the  region  R is  bounded  and  has  area  a,  then 

there  exists  an  absolute  constant  b such  that 

11m  iin  = b^a  almost  everywhere  (a.e.)  (2.2) 

n-»oD  yn  [ ] 

The  existence  of  the  limit  In  (2.2)  Is  part  of  the 
assertion.  The  constant  b has  been  estimated  by  Monte 
Carlo  experiments  to  be  .75.  Of  course,  the  rate  of 
convergence  of  the  sequence  L^/yn  Is  Important;  this 
aspect  of  the  problem  will  be  discussed  In  detail  In  Section 
4.3.  We  shall  also  discuss  there  the  case  In  which  the 
probability  distribution  Is  not  uniform. 

By  Theorem  1,  then,  the  length  of  the  shortest  path,  a 
random  variable.  Is  asymptotically  (i.e.  for  n large) 
equal  to  byaTn  with  probability  1.  The  value  bTas/n  Is  a 
; non-random  function  of  n and  so  for  large  n there  Is  no 

f distinction  between  the  random  and  non-random  versions  of 

I the  problem,  and  we  can  predict  with  probability  1 the 

(length  of  the  optimal  tour  through  any  random  points.  In 

particular,  suppose  that  an  algorithm  A yields  a tour 
A 

length  of  Ln  through  n random  points,  and  that  n Is 
very  large.  Then,  by  definition  (2.1),  this  algorltm  solves 
the  problem  to  within  Ln/bySyn  with  probability  1.  It  is 
In  this  asymptotic  probabilistic  sense  that  we  are  able  to 
investigate  the  optimality  of  various  algorithms. 


I 
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Conslder  the  following  algorithm  (similar  to  that  of 
Karp  [9])  which  yields  a path  for  any  Instance  of  the  TSP. 


Divide  the  region  R into  m subregions  each  of  area 
a/m,  and  label  them  r i , r2 , • . . rg, . Using  an  optimal 

algorithm,  construct  an  optimal  travelling  salesman  tour 
within  each  of  these  regions  individually.  Now,  visit  the 
regions  ri,r2,...rg,  in  order;  upon  completing  a tour  in 
region  r^,  visit  the  closest  unvlsited  point  in  region 
■"l  + l (take  rji^l  = r^)  and  thence  traverse  the  tour  in 
that  region.  [] 


For  n large  enough  there  will  be  n/m  points  in  each 

( v ) 

region;  the  tour  within  each  region  has  length  " (by 
Theorem  1) 


bln  I A » JiyETa 


(a.e.  ) . 


Calling  the  total  distance  Tq  we  have 


Tn  ^ mi^ynya  ■4- 


(a.e.  ) 


(3)  It  should  be  clear  that  these  equations  and  inequalities 
bold  In  the  asymptotic  probabilistic  sense  only.  The 
derivation  can  easily  be  made  rigorous.  The  appendix  of 
Chapter  II  will  illustrate  how  this  can  be  done. 
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(4 ) 

where  A is  the  diameter  of  the  region  R. 

Hence , 

Tn  ± bJa  + (a.e.) 

yn  Tn 

and,  given  any  c , for  n large  enough, 

^n  i bTa  + e (a.e.) 

7H 

So,  loosely  speaking.  Algorithm  1 is  'asymptotically 
optimal'.  More  formally,  we  have 

Corollary  J. 

Given  any  c>0,  there  exists  an  N(c)  such  that  for 
any  random  problem  of  size  n,  with  n2N(e),  Algorithm  1 
is  c-optlmal  with  probability  1.  [] 

It  can  further  be  easily  shown  - see  Karp  [9]  - that  if 
■ grows  as  loglogn  (so  that  m/n->0  as  n->oo)  and  an 
0(n.2'‘)  algorithm  (e.g.  that  in  [8J)  is  used  within  each 
subregion,  then  Algorithm  1 runs  in  time  O(nlogn)  a.e. 


(4)  The  diameter  of  R is  miniilx-yll;  x,yeR}.  The  term 
involving  the  diameter  becomes  negligible  in  the  limit.  The 
formulas  are  Independent  of  the  shape  (or  convexity)  of  the 
region.  For  example,  in  a region  with  shape 

(f" 

the  travelling  salesman  tour  will  still  be  hardly  likely  to 
exit  from  the  region  if  n is  large  enough. 
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Static . Single-Bus  £r.Qb.l.Stt8 

2.2.1  Optimal  Algorithms 

The  simplest  static  bus  problem  is  one  in  which  there 
is  only  a single  bus.  We  wish  to  choose  a feasible  tour 
through  the  n demand  pairs  so  as  to  minimize  the  total 
distance  travelled  (this  will  be  the  same  as  minimizing  the 
time  by  which  the  final  passenger  i.s  delivered).  This 
problem  differs  from  the  TSP  only  in  the  feasibility 
constraint . 

As  a first  observation,  note  that  the  problem  is 
NP-complete  (see  Observation  2.0  in  Appendix  II). 

(5 ) 

Let  us  define 

Lj]  s length  of  the  optimal  travelling  salesman  tour 
through  n random  points  in  R. 

Yn  s length  of  the  optimal  feasible  bus  tour  through 
n random  demand  pairs  in  R. 

Observe  that,  by  dropping  the  feasibility  constraint, 
we  would  obtain  a tour  of  length  L2n'  Hence 

L2n  i ^n  (2.3) 

An  upper  bound  to  Y^,  can  be  obtained  by  any  suboptimal 


algorithm  which 

might  be  suggested. 

For 

example,  if  we 

first  visit 

all  origins 

in  R 

and 

then  visit  all 

destinations,  we 

would  obtain 

a tour 

with 

1 2 

length  Ljj^Lj^, 

(5)  Both  Ln  and  Y^  are  random  variables. 
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1 2 

where  and  are  the  lengths  of  the  two  travelling 

salesman  tours  on  n random  points  each.  This  yields 

Yn  i 2Ln. 

We  can  easily  improve  on  this  upper  bound  as  follows. 
Divide  R into  two  equal  subregions  of  area  a/2  each, 
calling  them  regions  ri  and  r2 . First,  visit  all  origins 
in  r^  (there  will  be  n/2  points,  from  demands  of  the 


form 

(ri ,ri 

) and 

(ri,r2)  for 

n large).  Second,  visit 

the 

origins 

and  the  destinations 

from  demands 

(ri .r2) 

in 

»'2 

( there 

will  be 

3n/M  points 

here ) . Third  , 

visit 

*•1 

again,  visiting  the  remaining  points  there  - these  will  be 
the  n/2  destinations  from  demands  (ri,ri)  and  (r2,ri). 
Finally,  visit  the  last  destinations  in  r2  - there  wll'*  be 
n/«  points  from  demands  (r2*r2).  For  n large  this  tour 
has  approximate  length  (a.e.) 

b(yn/2  ♦ yTnTT  yn/2  )./a/2 

= 1.96byayn  * 1.96Ln. 

The  obvious  generalization  is  easy,  and  we  can  give  the 
following  algorithm,  producing  a feasible  tour  which  we  call 


Algorithm  2 (producing  the  tour  T^) 

Partition  the  region  R into  m equal  subregions, 
each  of  area  a/m,  and  label  them  ri,r2,...ra.  Visit  the 
regions  in  numerical  order;  in  each  region  rj^  collect  all 
the  origins  as  well  as  the  destinations  from  regions 
,r2 1 • • • >'i>  1 • This  is  the  'first  passage'  through  the 
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regions.  Following  this,  again  visit  the  regions 
ri,r2,...rj,  in  order,  visiting  all  the  remaining 
destinations  in  each.  This  is  the  'second  passage'  through 
the  regions.  Each  time  a region  is  visited,  the  tour  is  a 
travelling  salesman  tour.  [] 

It  is  shown  in  Lesma  2.2  that  the  tour  has  length 

jl.y?byayn  + yn0(1/m)  (2.4) 

3 

If  we  now  assume,  as  appears  reasonable  at  this  stage, 

that  there  exists  an  absolute  constant  c such  that 

lim  Xn  = c-Za  (a.e.)  (2.5) 

n->oo  Jn 

then  (2.4)  and  (2.5)  (noting  that  L2n  - -/Snbya  = -/SLu) 
together  give 

1 0 i 4y?b/3. 

More  difficult  is  Lemma  2.3  where  it  is  shown  that  the 
value  of  4y?b/3  is  also  a lower  bound  to  c (the  rigorous 
proofs  in  the  appendix  do  also  guarantee  the  existence  of 
the  limit  in  (2.5)  and  hence  too  that  c exists).  We  give 
a rough  outline  to  the  proof  of  Lemma  2.3>  Given  any 
optimal  tour  with  length  Yj),  it  is  possible  to  construct  a 
suitable  division  of  the  area  R into  m equal  subregions 
(for  any  m^l ) such  that  the  following  infeasible  tour  does 
yield  a lower  bound  to  the  value  Yn*  This  infeasible  tour 
visits  the  regions  successively,  twice  each,  in  exactly  the 
manner  of  the  tour  T^i  except  that  the  destinations  of  the 
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I 

form  (ri.rj),  i=1,2,...m,  are  visited  in  each  region  r^ 


on  the 

first 

passage 

through 

the 

regions 

and 

not  on  the 

second 

(6)^ 

Letting 

m-^CD 

the 

lengths 

of 

the  infeasible 

tours 

converge  to  4y7b./n/3, 

a . e 

From  Lemmas  2.2  and  2.3  we  can  state 

Ifafigrea  Z 

lim  Xil  = iiy?bya  a.e. 

n->oo  Vn  3 

= 1.89bya.  [] 

Henceforth  for  convenience  we  write 
I c = «y?b/3. 

He  shall  also  refer  to  the  'optimal  tour  T°'  with 

length  c/nVa,  where  T°  a lim  T^;  by  this  rough  statement 

m->OD 

we  mean  that  we  can  approach  arbitrarily  close  to  the 
optimal  value  of  by  a tour  given  by  Algorithm  2; 

^ this  is  still  asymptotic  in  n (with  n/m->OQ)  a.e. 

: More  precisely,  we  get  from  Lemma  1, 

i 


(6)  Note  that  if  there  are  n^  points  to  be  visited  in  a 
certain  region  r^^  and  if  n*  and  n”  are  visited  on  the 
first  and  second  passages  through  rj[^  respectively,  then 
the  tour  in  r^  has  length  behaving  as 

yj^  « bya/m(yirj  Vn^) 

with  ”i  * ®1  “ ®1 

The  value  yi  is  minimised  If  ln'-n"l  is  made  as  large  as 
possible.  ^ ^ e » 
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Glven  any  c>0,  there  exist  N(c)  and  m(c)  such 
that  for  n2N , the  tour  given  by  Algoritho  2 mlniaizes 


the  total  distance  travelled  to  within 
probability  1 . 


with 


2.2.2  CoBagpta  aJ^  Iht.  Xfiuc  !■ 


of  the  Tour  T 


Note  that  the  tour  T described  above  does  not 
uniquely  solve  the  minlaua  distance  problea.  Consider  a 
tour  Tg,  which  visits  the  regions  ( r i , r2  i • ■ • r^ ) In  the 
saae  aanner  as  on  the  first  passage  through  the 

regions.  On  the  second  passage,  however,  the  regions  are 
visited  In  the  order  (rB,rB.^ , . . .r^ } , and  once  again  all 
feasible  destinations  are  visited  In  each  region.  Clearly 
Tb  and  have  the  saae  length  (asyaptotlcally  a.e.) 


T''  > lla  Ti 

B^CO 

would  also  solve  the  problea. 


(7)  This  stateaant  can  be  strengthened  slightly:  Given  any 
OO  there  exists  an  H(t)  *nd  a function  N(b,c)  such 
that  for  all  >2H,  the  tour  alnlalses  the  distance  to 

within  (I't’C)  a.e.,  whenever  n2N(B,c). 
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In  later  sections,  particularly  for  the  dynanlc 

problem,  we  shall  be  interested  in  other  criteria;  there  is 

then  an  inefficiency  in  the  tour  as  it  is  'tail-heavy* 

with  passengers  having  on  the  average  a longer  travel-time. 

o 

This  Justifies  our  present  emphasis  on  the  tour  T^. 


There  is  yet  another  tour  which  minimizes  the  distance 
travelled.  Consider  visiting  regions  (ri,...rB,)  in  the 
same  order  on  two  separate  passages.  On  the  first,  collect 
in  each  region  r^  all  origins  of  the  form  C^'i.rj),  Jii . 
On  the  second  passage,  collect  all  origins  of  the  form 
(ri,rj),  J>i  and  simultaneously  deliver  all  feasible 
destinations  (these  will  be  all  of  the  destinations  in  rj^). 
It  is  easy  to  see  that  this  tour  will  have  the  same  length 
as  the  tour  T^,  but  will  be  even  more  tail-heavy  than  was 
T®’ 


We  might  emphasise  the  easy  principle  underlying  the 
optimal  tour  T^,  which  will  later  be  seen  to  be  important. 
Visit  the  regions  successively;  each  time  the  bus  exits  from 
a region  choose  as  next  region  the  one  with  the  most 
unvisited  feasible  points  (break  ties  arbitrarily),  and 
visit  all  of  these  points. 

aiMBla  iQura 

This  is  a suitable  Juncture  at  which  to  digress  to  make 
a further  important  observation  on  the  tour  T®.  It  belongs 
to  a class  of  tours  which  we  shall  term  'simple*. 
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Conceptually,  a aimple  tour  Is  one  that  visits  the  points  in 
R region-by-region.  The  tour  segment  in  r^  can  be 
constructed  at  the  time  at  which  the  bus  enters  ri,  in 
that  all  points  visited  in  rj^  were  feasible  at  this  time 
of  entry. 

More  precisely,  let  m be  any  Integer,  mil.  Define 
Sg,  the  set  of  'simple  tours  on  m subregions'  as 

follows.  A feasible  tour  T belongs  to  Sg,  if  there 

exists  a partition  of  R into  at  least  m subregions  of 
area  a/m  each,  with  the  two  properties: 

1)  if  T enters  a region  rj^  at  time  t^,  then  the  points 

which  are  visited  in  r^  were  already  feasible  at  time  t^^, 

2)  if  T visits  p points  in  r^,  then  we  can  assume  that 

these  p points  are  randomly  distributed  over  r^  (and  that 

T performs  an  optimal  travelling  salesman  path  on  these 
p points). 

Now  define 

CD 

S * U Si 
1*1 

Then  any  tour  in  S is  feasible,  and 

Si  c Si4,iC  S for  any  1*1,2,... 

We  call  S the  set  of  'simple  tours'. 

By  Corollary  2,  for  the  single-bus  static  problem  it  is 
enough  to  consider  only  simple  tours.  Most  tours  which  we 
encounter  will  be  simple;  in  Chapter  III  we  shall  restrict 
our  consideration  to  simple  tours  only. 
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2.2.3  Interesting  SubOBtlmal  AiKOrlthM 

It  is  of  interest  to  describe  briefly  some  suboptimal 
algorithms.  These  are  closely  related  to  currently  used 
techniques  and  will  be  needed  subsequently. 

The  following  might  be  considered  to  be  the  'Michigan 
Algorithm',  described  in  Chapter  I,  for  the  static 
single-bus  case. 


Partition  the 

region 

R 

into  m 

subregions . 

Calculate  an  optimal 

travelling  sa 

lesman  tour 

on  the  origins 

within  each  region. 

Visit 

all 

regions,  pe 

rformlng  this 

optimal  tour  in  each, 

and  linking 

the  regions 

together  using 

some  fixed-route,  of  any  length  g.  Then,  again  return  along 
the  same  fixed-route,  deviating  within  each  region  to 
deliver  all  destinations,  using  an  optimal  travelling 
salesman  tour  within  each  region. 

M 

Let  Yq  be  the  total  length  of  this  bus  tour. 

Then 

and 

v”  r- 

11m  in  s 2bya  a.e. 

n^oD  yn 

The  distance  g travelled  by  the  fixed-route  bus  is 
asymptotically  negligible.  It  will  henceforth  be  convenient 
to  think  of  the  region  R as  the  interior  of  a circle,  and 
the  subregions  as  sectors  (see  Figure  2.1). 
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Figur«  2.1  Partition  of  R into  subregions 

It  Is  then  easy  to  visualize,  for  exanple,  a bus  visiting 
the  subregions  in  any  required  order,  and  to  ellnlnate  the 
need  for  a flxed>route  bus.  However,  this  convenience  is 
not  necessary. 

The  'Michigan  tour*  described  above  is  now  seen  to  be 
only  a very  simple  extension  of  Algorltba  1 to  the  bus 
problem.  With  respect  to  distance  travelled,  it  can  be 
considered  the  same  as  the  tour  given  by  Algorithm  2. 

Most  Important,  from  the  practical  point  of  view,  is 
that  (see  (2.1))  the  algorithm  solves  the  problem  to  within 
2b/o  ~ 1.06  a.e. 

A second  suboptimal  algorithm  is  the  following 
' fixed>route'  algorithm.  Again  the  region  is  divided  into 
m subregions  and  a priori  it  is  specified  that  the  regions 
will  be  visited  first  in  the  order  r^,...ra 


end  then  in 
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the 

order  rni,r 

m- 1 » ' • 

•^1 

• 

On  the  first  passage,  origins 

of 

the 

form  (r^ 

UJ 

are 

collected 

and  destinations  of 

the 

form  ( rj , r 

j),  i< 

J 

ar 

e d 

elivered . 

In  region 

»*i. 

the 

re  will 

be 

(m- 

i + l 

) n/m 

2 origin 

s and  (1- 1 ) n 

/m2 

des 

tinations  vi 

sited , 

a 

to 

tal 

of  n/m 

points.  On 

the 

sec 

ond  passage 

all 

remal 

ning 

origins  a 

nd  destinations 

are 

vis 

ited  in  each 

regio 

n ( 

thi 

s wi 

11  be  fees 

ible) . 

F 

Let  Yq  be  the  total  length  of  this  bus  tour.  Then, 

Yn  = 2nibffi  a 
\m\m 

or,  more  precisely, 

F 

lim  Xa.  - 2bya  a . e . 

n->oo  Jn 

and  this  algorithm  also  solves  the  problem  to  within  1.0b 

a . e . 


2^  Static.  Multlple-Bua  Problems 

Introduction 


We  now 

consider 

the  static 

problem 

in  which  th 

ere 

are 

k buses. 

The  problem  rac 

ed 

is 

essentially 

one 

of 

multiprocess 

ing.  Any 

passenger 

can 

be 

served  by  any 

bus 

and 

we  are  to 

allocate 

the  passengers 

to 

buses  in  some 

optimal 

fashion. 
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Two  possible  generalizations  of  the  previous  criterion 
iomediately  arise,  and  we  get  the  problems: 

(i)  Minimize  the  total  distance  travelled  by  all  k buses 
(i.e.  minimize  fuel  used),  and 
(ii)  Minimize  the  time-to-completion  (i.e.  time-to-del 1 very 
of  the  final  passenger). 
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Lemma  3 . 1 


(8) 


lim  ^ a cVa 
n-^oD  ys 


a 


e . 


[] 


Thus,  by  increasing  the  number  of  buses  we  do  not 
improve  upon  the  total  distance  travelled  or  upon  the  fuel 
consumption.  If  instead,  the  final  time-to-delivery  of  all 
passengers  is  of  interest,  we  must  Investigate  Y|,.  Then  we 
require  that  each  of  the  buses  absorbs  part  of  the  load,  and 
the  Jobs  of  collecting  and  delivering  passengers  must  be 
executed  in  parallel.  For  efficiency,  all  passengers  should 
travel  an  equal  distance,  so  that  all  buses  are  busy  all  the 
time.  Clearly  too,  we  must  have 


Y 


k 

n 


2 


k 


We  ask,  under  what  circumstances  do  we  in  fact  get 


k 


Theorem  3 to  follow  indicates  that  this  lower  bound  is 

( 9 ) 

attainable  (in  the  limit)  if  passengers  may  transfer 
between  buses,  as  in  the  scheme  at  Michigan.  If  passengers 
cannot  transfer,  so  that  the  bus  collecting  a passenger  must 
also  deliver  him,  as  in  the  scheme  at  Rochester,  then  each 


(8)  The  proof  of  Lemma  3.1  which  is  given  in  the  appendix  is 
somewhat  Indirect,  and  an  easier  proof  has  been  found  to  be 
evasive.  However,  the  proof  does  Identify  additional 
problems  as  being  of  interest:  these  are  problems  in  which 
at  moat  i transfers  are  allowed  during  the  tours  of  the 
buses,  for  i«1,2,3,... 

(9)  At  a 'transfer  point',  at  least  two  buses  meet  and  the 
passengers  on  board  are  able  to  move  arbitrarily  between 
them. 
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bus  must  travel  a greater  distance. 


2.3.2 


Theorem  ^ 


Allowed 


If  transfers  are  allowed 

yk 

lim  in  = fiya 
n-^oD  yE  k 


Consider  the  following  algorithm,  producing  a tour 
k° 

which  we  call  Tg  . 


I 

i 


Figure  2.2  Algorithm  3 


Algorithm  2 (producing  Tg  ) 

Let  n be  any  positive  integer.  Partition  the  region 
R Into  mk  subregions  each  of  area  a/km.  Bus  1 will 
be  visiting  only  regions  >^(1>1  )a4-1  • • • 'rimt  l>1,...k. 

First,  let  each  bus  1 visit  region  r ( , collecting 
all  origins  there.  Then,  let  all  buses  meet  at  a 
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prearranged  transfer  point  (the  point  0 in  Figure  2.2)  and 
passengers  transfer  onto  the  bus  that  will  ultimately  be 
visiting  the  region  containing  the  required  aestination. 
Each  bus  i then  visits  region  >”(i-l)m+2>  collecting  all 
origins  tnere  and  delivering  all  feasible  passengers  (from 
r(j-l)m+1»  J = 1. .••!£).  The  k buses  meet  at  the  next 
transfer  point  when  passengers  who  have  Just  been  collected 
again  transfer  to  their  required  buses.  In  this  way,  after 
m transfers,  the  buses  will  have  collected  all  origins. 
Then  (in  the  manner  of  Algorithm  2)  bus  i revisits  regions 
'"(i-l  )m+1  I ■ • •*‘im»  delivering  the  remaining  passengers.  [] 

The  m transfers  enable  a passenger  to  be  on  the  bus 
visiting  his  aestination  as  soon  as  possible.  Whenever  a 
bus  enters  a region  at  time  t,  it  is  able  to  visit  every 
destination  tnere  that  was  feasible  at  t,  irrespective  of 
which  bus  visited  the  corresponding  origin.  Each  bus  visits 
2n/k  points  and  effectively  serves  n/k  passengers. 
Furthermore,  each  bus  travels  through  (1/k)th  of  the  total 
area  a (also,  each  passenger  needs  to  travel  at  most 
through  an  area  a/k)  and  travels  a total  distance 
cytiTTcyiTTc  » cyn/a/k.  The  following  corollary  is  then 
obvious  (but  see  the  detailed  proof  of  Theorem  3). 

CarailarY  1 

Given  any  c>0  tnere  exist  N(c)  and  m(£)  such  that 

IrO 

for  n2N  the  tour  T^  given  by  Algorithm  3 is  c -optimal 
for  problem  (11)  (i.e.  minimizing  completion-time)  with 
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probability  1 . [ ] 

For  the  tours  given  by  Algorithm  2,  we  defined  a 

'tour'  T°  = 1 im  There  was  no  conceptual  difficulty  in 

m-^oD 

doing  this:  the  larger  m became  the  sooner  the 

destination  whose  origin  had  Just  been  collected  became 

k® 

feasible.  The  tour  Tn  of  Algorithm  3 now  requires  m 
transfer  points;  in  the  analogous  limit  as  m->co,  a limiting 
optimal  tour  would  require  continuous  transfers  on  the  first 
passage  through  the  region  R.  We  thus  refrain  from 

defining  a corresponding  limiting  'tour'. 

If  we  define  the  'Michigan  tour  with  k buses'  to  be 

k© 

the  tour  yielding  T^  (so  this  tour  has  only  one  transfer) 

what  do  we  lose?  As  before,  the  distance  travelled  by  each 

bus  before  the  transfer  (i.e.  n/k  origins  in  an  area 

a/k  are  collected)  is  bs/n/kVa/k . After  the  transfer, 

n/k  destinations  are  visited  in  the  same  area.  Hence,  if 
k” 

is  the  completion  time  which  this  algorithm  yields, 

then 

kM 

lim  ^n  s Ziiy/a  a.e. 

n->0D  yE  k 

Theorem  3 gives  that  this  tour  solves  problem  (li)  to  within 
(2b/k)/(c/k)  s 1.06  a.e.  once  again;  we  lose  but  6^  by 
restricting  ourselves  to  only  one  transfer. 
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It  is  possible  to  improve  on  this  algorithm  slightly. 
Passengers  with  both  origin  and  destination  in  the  same  area 
serviced  by  bus  1,  say,  might  in  fact  be  able  to  be 
delivered  before  the  transfer  point.  Thus,  divide  the 
region  into  mk  subregions,  with  each  bus  servicing  m 
regions,  and  a single  transfer  talcing  place  after  the  first 
passage  through  these  regions.  Now  m can  be  made  large 
with  no  conceptual  di f f iculty ^ ^ ^ ^ . Such  schemes  will  be 
discussed  in  greater  detail  for  the  dynamic  problem  of 
Section  3.3. 

^3.3. 3 ILa  Iranalera  All0Mfi.d 

We  now  turn  our  attention  to  schemes  in  which  transfers 
are  not  allowed.  Of  course.  Lemma  3.1  still  nolds,  but 
Theorem  3 does  not.  It  is  easy  to  see  that 

Tn  2 li  (2.6) 

k 

and,  for  k>1,  it  appears  that  a strict  inequality  holds  if 
there  are  to  be  no  transfers.  We  do  not  have  results  as 
strong  as  Theorem  3 for  this  case. 


(10)  This  scheme  yields  a completion-time,  satisfying 

lim  iSH  « lbv^[(1  ♦ 1/k)^''^  - (1  - a.e. 

n^oo  yn  3 

The  proof  of  this  is  similar  to  that  of  Theorem  6 (Appendix 

1X1). 
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( 1 ) Two-Bua  Caae 

For  k=2  consider  the  following  algorithm,  which  is  a 
'fixed-route  scheme'  for  two  buses  (recall  the  fixed-route 
algorithm  with  one  bus).  Once  more  divide  R Into  m 
subregions  ri,r2,...rD'  bus  I visit  the  regions  In 

order  ( r i , r2 , • . . rn ) collecting  and  delivering  all 

passengers  with  demands  ^^lirj),  J>i.  Let  bus  II  visit 
the  regions  in  order  ( rn , rj,.  i , . . . r i ) , collecting  and 
delivering  passengers  with  demands 

Passengers  with  demands  (ri,ri)  are  divided  evenly  between 
the  buses  and,  as  we  can  let  m->ao  these  passengers  will  not 
cause  difficulties.  (See  Figure  2.3.) 


Figure  2-3  Two-bus  static  scheme  without  transfers 

Each  bus  will  visit  n/m  points  in  each  region  as 
before:  for  example,  bus  I in  region  r^  will  collect 

I 

2 2 
(■-1)n/m  origins  (ri,rj),  j>1,  n/20  origins 
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2 

(ri,ri)  and  will  deliver  n/2m  destinations  (ri,ri). 

the  completion-time  which  this 

mb  IT  a 
jmjm 

by¥  a . e . 

Having  considered  numerous  other  two-bus  schemes,  we 

are  lead  to  conjecture  that 

lim  Xll  = bTa  a.e. 

n->aD  yn 

i.e.  that  this  scheme  is  asymptotically  optimal  (a.e.)  for 
the  minimal  completion-time  problem  (problem  (ii))  with  no 
transfers . 

Be  this  as  it  may,  from  (2.6)  we  obtain 

v2  „1 

lim  iQ  2 lim  ^ n 
n->CD  ya  n->OD  2yn 

s fiTa  a.e. 

2 

so  that  the  fixed-route  scheme  at  worst  solves  problem  (ii) 
with  no  transfers  to  within  b/(o/2)  s 1.0b  a.e. 

l ii ) Three-Bus  Case 

A similar  scheme  which  we  conjecture  to  be  optimal  for 
k>3  buses  is  ' ^e  following.  (See  Figure  2.4.) 


2^ 

Thus,  if  Yn 
algorithm  yields,  then 

„2F 

lim  ln_  = 
n->0D  yn 
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Rl 


R2 

Figure  2.4  Three -bus  static  scheme  without  tronsfers 


Divide  R into 

3a  subregions 

PI . 

P2.' 

* • ^3b  * 

m 

2m 

3m 

Rl  3 U 

Pi.  ^2  * 

0 ri, 

, R3  = 

U 

Pi- 

Tnen  bus 

1=1 

lsm-«-1 

i = 2m4'1 

services 

demands 

(Ri.Rl), 

, (Rl. 

R2) 

and 

(R2.R1) 

travelling  in  order  through  ( r i , r2 , . . • r^ ) , through 

(r'2B  f »'2m- 1 » • • • *'m+1  ^ and  through  (r  i , r2 , . . . rn)  again. 
(Call  these  the  first,  second  and  third  passages 
respectively.)  On  the  first  passage  the  bus  collects  in 
Pi.  IlliBi,  all  origins  of  the  form  (Pi.rj).  iiJi2m,  and 
delivers  all  feasible  destinations,  l.e.  of  the  form 


(rj,rj),  i.lJ<i.  On  the  second  passage  in  regions 
Pi»  m+1iii2m  (in  reverse  order!)  origins  of  the  form 
(ri,rj),  ii.J<i  and  all  feasible  destinations 

(Pj.Pi).  tlJlB  and  i<Jl2m  are  visited.  Finally,  on  the 
third  passage,  all  remaining  destinations  in  r^^,  tll^m  are 
visited,  i.e.  (Pj.Pi)i  ■•»1U.i2m 


and 


(Pl .Pi) . 


At  the 


f 
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saoe  time,  bus  II  services  demands  (R2>^2)i  (R2>^3) 

and  (R3,R2)i  and  bus  III  services  demands  (R3,R3), 

(R3,Ri)  and  (Ri,R3),  in  a similar  way. 

As  m-^co,  the  distance  travelled  by  each  bus,  defined 


to  be  Y: 


satisfies 
3^ 


lim 

n->ao  yn 


a . e . 


. 798bya 

We  derive  this  fact  in  the  appendix,  in  Proposition  3.2. 


mil  It-Bua  JCflSS 


! The  approach  taken  above  can  be  used  to  obtain  what 

I 

appear  to  be  good  schemes  for  k buses,  with  k large. 


2 

Suppose  there  are  kss  buses.  Divide  R into  s 

subregions  and  let  each  bus  service  a particular 

demand-type,  i.e.  i i 1 . 2 , . . . a . For  iWJ,  the 

2 

bus  serving  (ri,rj)  will  collect  n/s  points  in  each  of 

the  regions  rj^  and  rj,  yielding  a time-to-dellvery  of 

2b  _n_  a s 2bJriVB  (asymptotically  a.e.)  (2.7) 

• l<3/« 

The  buses  serving  demands  (rj^.r^}  would  complete  their 
tours  in  less  time,  so  (2.7)  is  the  maximum 
time-to-dellvery. 


If  kBs(s-l),  we  are  able  to  do  better  since 
then  possible  to  have  all  buses  busy  all  of  the  time 
with  s subregions  in  R,  let  each  bus  service  a 
IWj  and  1/2(s-1)  of  each  of  the  demands 


it  is 
Again 
demand 

^*'i.»*l) 
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and  (ri,rj).  Then  In  each  of  its  regions,  a bus  will  visit 

1 JL  _IL  = n points. 

(s-1)  .2  2 s(s-1) 

S 9 

The  total  distance  travelled  by  each  bus  is  then 


2b  a n 
J sj  s ( s-1  ) 

2na 


t1/2 


Lk(yTnr7-1)J 


(a . a . e . } 


(2.8) 


For  very  large  k the  values  in  (2.7)  and  (2.8)  are 
not  of  too  much  interest.  However,  it  is  believed  that  the 
values  are  then  'almost  optimal':  they  indicate  the  minimal 
tlme-to-delivery  as  a function  of  k,  for  k large. 


The  schemes  for  the  multi-bus  problem  that  we  have 
analysed  are  seen  in  retrospect  to  have  a certain 
specialization  of  the  buses  in  common.  For  the  case  in 
which  transfers  are  allowed,  the  buses  are  restricted  to 
specialist  regions . For  the  case  in  which  transfers  are  not 
allowed,  the  buses  specialize  in  collecting  demands  of  a 


certain  type  only 


CHAPTER  111 


Iil£  AVERAGE  FLOW-TIME  CRITEBIQW  AHI)  IM  MNAHIC  £A££ 

Two  Inadequaclea  of  the  models  In  Chapter  II  are  the 
distance  criterion,  which  is  not  really  relevant  to  an 
individual  passenger,  and  the  static  nature  as,  in 
actuality,  we  must  consider  an  ongoing  process,  with  demands 
arriving  over  time. 

This  chapter  investigates  these  aspects  of  the  problem 
by  generalizing  the  previous  results.  In  the  dynamic  case 
we  shall  consider  the  average  flow-time  criterion,  so 
Section  3.1  investigates  this  criterion  for  the  static  case. 
Then,  Section  3.2  analyses  the  (single-bus)  dynamic  problem 
and  its  relationship  with  the  static  problem.  Finally, 
Sections  3.3  and  3>4  study  the  multiple-bus  dynamic  problem. 

Once  again,  details  of  proofs  are  collected  in  an 
appendix.  The  notation  and  assumptions  are  the  same  as 
those  described  in  Section  2.0. 

1 Iil£  Static.  Sinale-Bus.  Averaae  Flow-Time  Problem 

r 

The  criterion  used  in  Chapter  II  - that  of  minimizing 
the  t Ime-to-delivery  of  the  final  passenger  - is  not 
necessarily  satisfactory  for  an  individual  passenger.  For, 
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he  is  Interested  rather  In  his  personal  wait-time  and 
travel-time.  Therefore,  consider  minimizing  the  average 
flow-time,  i.e.  the  average  t ime-to-delivery , of  all 
passengers.  This  criterion  will  be  used  in  the  dynamic 
problem;  the  present  section  is  needed  as  a link  between  the 
static  and  dynamic  cases. 


We  need  some  additional  notation.  For  any  given 
problem  instance  of  size  n,  let 
Oq  - set  of  n origins 

= set  of  n destinations. 

Any  feasible  tour  T defines 

t(oilT)  s time  at  which  origin  o^  is  visited 
tCd^lT)  s time  at  which  destination  d^  is  visited. 

Let 

n 

Wn(T)  » 1 I t(oilT)  s average  waiting  time  of  T 
n i*  1 


fn(T) 


n 

J.  £ t(dj.lT)  z average  flow-time  for  all 
n izi 

destinations  of  T. 


The  problem  is,  then,  given  a set  of  n demand  pairs 
distributed  at  random  in  the  region  R with  area  a, 
minimize  f,j(T),  with  respect  to  all  feasible  tours. 

We  shall  be  solving  a simpler  version  of  this  problem: 
essentially,  we  discretize  the  tours  by  considering  only 
simple  tours  in  S,  those  that  visit  the  points  in  R 
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reglon-by-reglon^ ^ ^ . Then,  decisions  are  made  at  discrete 
Instants  of  time  - whenever  a bus  exits  from  a region.  The 
choice  of  a tour  corresponds  to  the  choice  of  a strategy:  at 
each  stage  that  a bus  exits  from  a region  we  must  declde:- 
(1)  which  region  is  to  be  visited  next, 

(11)  which  points  are  to  be  visited  in  that  region. 

The  restricted  problem  is 

minimize  {fn(T),  TCS)  (3.1) 

In  order  to  provide  a Beardwood-type  result  for  problem 
(3.1)  we  first  consider  a problem  with  an  easier  solution: 

minimize  {Wn(T)  ♦ fn(T),  TCSn)  (3-2) 

We  draw  on  classical  scheduling  theory  for  our  Informal 
discussion.  Becall  the  following  result  (see  e.g.  Conway 
[6]).  Given  a single  machine  and  n Jobs  to  be  processed, 
the  average  flow-time  is  minimized  if  the  Jobs  are  ordered 
in  increasing  processing-time  order. 


Suppose  that  each  of  the  m subregions  has  area  a/m, 
and  that  nj,  points  are  visited  in  a certain  region.  Then, 


the  time  spent  there  is,  for  large  ni i approximately 


(1)  Recall  the  definition  in  Section  2.2.2. 

It  might  aeem  that  the  restriction  to  S is  a stringent 
one.  However,  for  the  ainiaum  distance  problems  in  Chapter 
II  there  is  no  loss  of  generality  if  we  restrict  ourselves 
to  S,  as  we  are  interested  in  c -optimal  solutions:  for 
any  optimal  tour  T*  there  exist  an  H and  a T,^  in  S,. 
aueh  that  " " 

T(Th)  - I(T*)  < c a. a. 

for  n large  enough.  It  is  believed  that  this  observation 
is  valid  for  the  average  flow-time  problem  too,  but  a 
rigorous  proof  proaisea  to  be  awkward. 


i 
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bynX^a/m , and  we  can  take  the  time  apent  on  each  of  the 
points  to  be  bTa/ymynX.  (Recall  that  the  n^  points  are 
randomly  distributed  In  the  region  and  nj^  Is  large,  so  we 
can  assume  that  the  time  spent  between  visitations  Is 
constant.)  Now,  to  minimize  the  flow-time  of  all  points . the 
scheduling  result  described  above  Implies  that  we  must  visit 
first  the  points  which  will  require  least  time,  l.e.  we 
must  minimize  bya/yEynJ.  To  do  this,  use  the  following 
rule . 

AlgflritbB  1 

Partition  the  region  R Into  m subregions  of  equal 
area.  At  each  stage,  upon  exiting  from  a region,  choose  as 
next  region  that  with  the  most  feasible  points,  and  visit 
all  feasible  points  there.  [] 

This  procedure  corresponds  precisely  with  the  tour  T^ 
given  by  Algorithm  2 (ties  are  broken  arbitrarily).  It 
minimizes  the  average  flow-time  of  all  points . l.e.  origins 
and  destinations,  which  Is  the  criterion  of  problem  (3*2). 

Observe  that  the  tour  T°  has  a long  average  waiting 
time,  Wn(T^):  passengers  wait  for  collection  while  as  many 
previous  passengers  as  possible  are  delivered.  Since  T° 
minimizes  Wn(T)  * fn^T)  in  (3.2),  we  might  expect  that 

works  well  at  minimizing  ^n^^)  in  problem  (3*1). 
This  Is  Indeed  the  case,  and  T°  Is  asymptotically  optimal 
in  (3>1}  again.  This  Is  stated  and  proved  (In  the  appendix) 


In  Theorem  4. 
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! 

i 


I 


> 


More  precisely,  in  Leoma  4.1  we  aerlve  the  facts  that, 
for  the  tour  given  by  Algorithm  2, 


lim 

Wn ( Tn  ) - 

.575bya  ♦ on/m) 

a.e. 

n->OD 

yn 

lim 

f n ( Tu  ) - 

1 . I42bya  + 0( 1/m) 

a.e. 

n->a> 

yn 

Hence  we  have,  from  the  considerations  above, 

Leami 

Let  Wn  fn  = lnf{Wn(T;  ♦ fn(T),  TCS} 

Then , 

lim  Wn  fn  s a.e. 

n->oo  yn 

Further,  given  any  c>0,  there  exists  an  M(c)  such  that 

for  m^M  and  n large  enough,  Tg  is  c-optimal  for 

problem  13.2)  a.e.^^^  [] 

IbBorea  il 

Let  fj  = inf{fn(T),  TCS). 

Then, 

lim  In  * 1 . a.e. 

rv*co  -/a 

Further,  given  any  c>0,  there  exists  an  M(c)  such  that 

for  m2.M  and  n large  enough,  ia  c-optimal  for 

problem  (3.iJa.e.  [] 

Thus,  we  have  that  the  tour  T°  la  nicely  robust.  It 
alnlaixes  a oost  relevant  both  to  the  passengers  and  to  the 
operator  of  the  system. 

(2)  Reoall  the  footnote  at  Corollary  2.  We  require  this  ) 

stronger  statement  for  the  proof  of  Theorem  4. 
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It  is  worth  observing  that  by  combining  Lemma  4.2  and 
Theorem  4,  we  obtain  the  result  that  the  tour  T°  minimizes 
with  respect  to  Tes  any  criterion 

g.Wn(T)+fn(T),  g 6 [0,1]. 
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^.2  The  Pvnamlc  Case ; Single-Bus  frO^l.CM 

Possibly  the  most  glaring  inadequacy  of  the  previous 
problem  formulations  was  the  assumption  that  all  demands 
were  known  at  the  starting  time,  t=0.  For  a more  realistic 
formulation  we  must  assume  that  the  demands  arise  as  time 
proceeds  according  to  some  Intermittent  arrival  process.  In 
this  section  we  study  the  extension  of  the  results  from  the 
static  case  to  this  dynamic  version. 

As  before,  we  shall  assume  that  any  particular  demand 
has  both  origin  and  destination  drawn  independently  from  a 
uniform  probability  distribution  over  the  region  R and 
that  the  demands  are  all  Independent  of  one  another.  For 
simplicity,  we  assume^^^  that  demands  arise  at  a constant 
rate,  q. 

To  complete  the  problem  specification  we  need  a 
performance  criterion.  It  is  clear  that,  to  minimize 
distance  only  without  regard  to  time,  is  uninteresting:  we 

can  always  wait  until  next  year  and  tnen  collect  and  deliver 
all  demands  together.  So  we  choose  to  minimize  the  average 
flow-time.  Thus,  if  a particular  demand  arises  at  time  t^ 
and  the  destination  is  visited  at  time  t^*  the  flow-time 
for  this  demand  is  is  averaged  over  all 

demands. 

(3)  This  assumption  is  unnecessarily  harsh,  and  certainly  we 
could  assume  a Poisson  arrival  process  with  mean  q. 
However,  the  assumption  is  consistent  with  our  informal 
presentation  and  it  eliminates  the  need  for  many 
'expectations' . 
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As  In  Section  3>1  we  discretize  the  problem  by 

restricting  ourselves  to  the  class  of  simple  tours,  S. 

A standard  easy  heuristic  for  solving  the  dynamic 
version  of  a decision  problem  Is  to  use  the  solution  of  the 
static  version  recursively  - this  loosely  corresponds  to  the 
'open-loop-optimal  feedback'  control  scheme  described  by 

Dreyfus  [7].  For  our  problem  the  technique  Is  the 

following.  At  each  stage  compute  the  optimal  (minimal 
flow-time)  static  tour  for  visiting  all  remaining  points, 
assuming  that  no  new  demands  will  arise.  Use  this  tour  to 

yield  the  Initial  decision  for  the  dynamic  problem  - l.e. 

which  region  to  visit  next,  and  which  points  there.  By  the 
time  this  has  been  Implemented,  new  demands  will  have 

arisen.  With  this  new  Initial  state,  recompute  an  optimal 

static  solution. 

In  general,  the  heuristic  Is  suboptlmal  as  It  does  not 

i 

take  Into  account  that  future  demands  will  arise.  But, 

consider  Its  application  to  our  problem.  We  claim  that  the 
static  minimum  flow-time  problem  which  Is  faced  at  any 
particular  stage  Is  solved  by  the  same  principle  given  In 
Algorithm  4:  at  each  stage  visit  the  region  with  most 

feasible  points  and  visit  all  such  feasible  points  there. 
We  do  not  prove  this  claim;  a detailed  proof  would  be  as 
long  as,  and  more  Involved  than,  the  previous  proofs,  and  no 
additional  Insights  would  be  gained. 
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The  following  is  a translation  of  the  heuristic  into  a 
proposed  solution  for  the  dynamic  problem  - Algorithm  4 is 
applied  recursively.  We  call  the  tour  it  is  believed 
that  no  confusion  with  the  static  tour  can  arise. 


Algorithm  5.  (yielding  the  tour  T^) 


Partition  the  region  R 

into 

m 

subregions  of 

equal 

area.  In 

each  subregion 

which 

is 

visited,  perform  a 

travelling 

salesman  tour  on 

all 

feasible  points. 

Upon 

exiting  from  a subregion,  choose  as  the  next  one  that  with 
most  feasible  points.  [] 

In  order  that  the  proposed  static  tour  at  each  stage  be 
(almost)  optimal,  the  arrival  rate  q must  be  very  large. 
Then  the  tour  will  visit  the  regions  consecutively. 
For,  with  probability  1,  the  region  with  the  most  feasible 
points  will  also  be  the  one  with  longest  elapsed  time  since 
it  was  last  visited. 

How  well  does  Algorithm  5 perform?  Using  a fundamental 
result  from  queuing  theory  we  argue  the  following  lemma  in 
the  appendix. 

Lemma  5 . 1 


If  the  system  is  in  steady-state,  then  with  probability 
1,  the  tour  given  by  Algorithm  5 minimizes,  among  tours 
in  Sg,  the  average  flow-time  of  all  passengers.  [] 
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Wlth  this  knowledge  it  is  possible  to  give  a 
Beardwood-type  relationship  for  the  dynamic  problem  when  it 
has  reached  an  equilibrium.  Since  q is  the  arrival  rate 
of  passengers,  in  order  that  the  system  be  in  a stable 
steady-state,  we  must  serve  q passengers  per  unit  time. 
Thus,  on  average,  2q  points  must  be  visited  per  unit  time. 
For  the  tour  in  Sn,  described  above,  let 

dq  s the  time  spent  in  each  subregion  when  the  arrival 
rate  is  q. 

During  this  time,  2q6q  points  must  be  visited;  so  by 

Theorem  1 (recall  that  the  buses  travel  at  unit  speed), 

lim  r bl^  a.e. 

q->ao  y^qSo  m 


lim  ^ = 2b^a 

q->GO  q m 


(3.3) 


Define  P-  s mSg  as  the  'period* 


this  is  the  time 


required  for  the  bus  to  perform  a circuit  around  R, 


visiting  all  subregions.  Then, 

lim  Is  X 2b^a 
q->0D  q 


(3.'*) 


From  this  it  is  easy  to  derive  the  formula  below  (see 


the  appendix). 


IfaBorm  1 


Let  Fq,m  optimal  average  flow-time  for  tours 

In  Sgi  when  the  arrival  rate  of  passengers  is  q.  Then, 


lim  X 2b^i(U1/m) 

q*>aD  q 


(3.5) 
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► 


Hence,  if  Fq  Is  the  optimal  average  flow-time  for 
tours  in  S,  then, 

lim  Is  = 2b^a  a.e.  (3-6) 

q->OD  q t ] 

This  theorem  is  analogous  to  Theorem  4,  but  for  the 
dynamic  problem.  It  is  trivial  that,  by  restricting 
ourselves  to  Sg)  instead  of  to  S,  we  solve  the  minimum 
flow-time  problem  to  within  (1-*-1/m)  a.e. 


3 Dynamic 


PXAtlleBl 


allowed 


Now  let  us  extend  the  last  section  to  the  k-bus 
problem . 


To  Introduce  the  approach,  consider  a 'Dynamic  Michigan 

k° 

Scheme',  a tour  which  we  call  (again)  T^  . 

ilBOrlthB  .LJ.  (yielding  t![°) 

Partition  the  region  R into  k subregions,  each  of 

( 4 ) 

area  a/k.  Let  all  buses  meet  at  a common  transfer  point 
at  times  6,20,30,.<..  During  time  lie,(i-f1)0]  each  bus 
visits  its  own  subregion  collecting  passengers  who  arrived 
during  [(1-1)6,01  and  delivering  passengers  who  arrived 
during  [ (i-2 )0, ( i-1 )e ] . At  the  transfer  point  each 
passenger  transfers  onto  the  bus  serving  his  required 
destination.  [] 

(A)  Recall  Section  2.3.2  where  this  was  justified  by  Figure 
2.2.  We  omit  any  mention  here  of  a 'line-haul'  bus. 
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Here,  6 s dq,  the  steady-state  time  spent  by  each  bus 
in  its  region.  During  this  time  q6q  additional  demands 
arise  with  2qeq/k  new  points  in  each  region.  Then,  as 
before , 


1 . e . 


11m  s b a 

q*»aD  y2qeq /k  k 

lim  ^ = 2b^a 

q-^OD  q ^^2 


The  period  Pq  equals  8q , so 

11m  Is  = 2b^4 

q-^CD  q j^2 


a . e . 


a . e . 


a . e . 


Now  seek  the  average  tlme-to-delivery . Consider  a 

passenger  arriving  during  time  [O.Oq]  where  0 and  Sq 

are  transfer  times.  This  passenger  is  collected  during 

[eq,28q]  and  is  delivered  during  [28q,38q]«  The  average 

time-of-arrival  is  8q/2,  the  pick-up  time  is  on  average 

38q/2  and  t>''  delivery-time  is  on  average  58q/2  (the 

expected  wait-time  is  8q,  as  is  the  expected  travel-time ) . 

M 

So  the  average  flow-time  by  this  algorithm,  viz.  ^q,ki 
28q  satisfying 

pM 

11m  « 4b*^a  a.e.  (3.7) 

q->OD  q ^2 

With  ks1  this  tour  is  precisely  T°  and  thus  (3.7)  is  but 
equation  (3.5)  with  msl,  as  it  should  be. 

Ic® 

The  tour  T^  is  in  A demand  arising  during 

[0,8q]  must  wait  until  the  next  period  before  it  can  be 
collected.  This  wait-time  can  be  Improved  by  further 
dividing  the  region  of  eech  bus  (with  area  a/k)  into  m 


f 

4 
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equal  subregiona  of  area  a/km  each.  Let  each  bus  visit 

all  of  its  m subregions  sucessively,  visiting  as  many 

origins  as  possible  in  each,  and  all  the  feasible 

destinations  - now  'feasible*  refers  to  the  fact  that  the 

corresponding  origin  has  oeen  visited  gjxsl  the  passenger  is 

on  the  bus.  Alter  visiting  tneir  m subregions,  the  buses 

all  meet  at  the  same  transfer  point  as  before,  passengers 

transfer,  and  the  cycle  is  repeated.  We  call  the  resulting 
)( 

tour  it  is  in  Sm);  and  it  requires  a single 

k® 

transfer  per  period.  Compared  with  , the  waiting-time 

of  passengers  is  reduced  and  passengers  with  origin  and 
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But,  as  pointed  out  earlier,  it  may  be  of  interest  to 
limit  the  number  of  transfers  per  period.  Suppose  we  insist 
on  only  one  transfer  per  period,  and  consiaer  tours  in  S|g^. 
Then  T , is  optimal  I 

1 I B 

More  generally,  consiaer  the  k-bus  problem  in  which  we 
have  at  most  h transfers  per  period,  and  we  allow  tours  in 

Skhm*  Then  we  have  the  following  algorithm  yielding  the 

. _k 

tour  Th,iB- 

Aitorltbl  (producing  Th^o,} 

Divide  the  region  R into  knm  subregions  of  equal 
area.  Designate  hm  subregions  for  each  of  the  k buses 
and  label  them  **i  * **2  * ' ’ ’ ''hm  ’ Let  bus  i visit 

successively  its  regions  ( r^  , r^  , . . . rf; ) . Following  this, 
all  buses  meet  at  a preassigned  transfer  point.  Tnen  bus 

i visits  regions  ( r^  and  a second  transfer 

occurs.  After  h such  transfers,  the  buses  again  visit 
regions  and  the  cycle  repeats. 

Every  time  a bus  visits  a region  it  visits  as  many 
points  as  possible  in  that  region;  these  are  all  the 
uncollected  origins  and  all  the  feasible  destinations  of 
passengers  already  on  board.  [] 

liLtoraa  £ 

The  tour  Th,m  optimal  (asymptotically  a.e.)  for 
the  k-bua  minimal  flow-time  problem  when  tours  are 
restricted  to  and  only  h transfers  are  allowed  per 
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period . 


Let  Fq;h,ni  t^®  average  flow-time  resulting  from 
tnis  tour.  Then, 


pk 

lia  *^<j;h,m  . 

_L(  Ifl/k) 

♦ Kl-l/k)! 

a.e. 

(3.9) 

q->OD  q 

“ L 

mb 

h J 

[] 

From  (3. 

b),  when  there  is  a 

single  bus, 

the 

optimal 

average  flow- 

2 

time  is  2b  aq.  From 

13.9),  When 

there 

are  k 

buses,  the 

conceptual 

optimal 

tour  in 

S = Soo 

with 

continuous 

transfers 

(l.e.  h->ad  has  a 

flow-time  of 

2 2 2 

2b  aq/k  , and  we  have  a k -fold  Improvement.  This  is 

(5) 

analogous  to  the  static  case  ; each  of  the  k buses 
visits  (1/k)th  of  the  total  area  a and  has  an  effective 
arrival  rate  of  q/k. 

Using  Equation  (3.9J  it  is  possible  to  investigate  the 
optimality  of  suggested  algorithms  with  respect  to  any 
restricted  problem  formulation. 

For  example,  among  tours  in  which  only  a single 
transfer  is  allowed  per  period,  what  is  lost  by  using  the 
tour  T®  7 ¥all, 

fq(T®)  « a.e. 

Letting  m^oD  with  ba1  we  get  from  (3.9)  tffat  the  optimal 
tour  oas  a flow-time  of  b aq(3-1/k)/k  and  so  solves 

the  problem  to  within 

(5)  Zn  the  atatlo  case  In  Section  2.3.2,  the  total  distance 
travelled  was  (asymptotically)  proportional  to  ^nvii, 
yielding  a k-fold  improvement. 
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[‘i/3,2] 


- 1 ♦ 1 e 

3-1/k  3-1/k 

and  we  lose  at  least  33^{.  Tnls  is  considerably  more  than 
the  loss  of  6}  obtained  for  the  static  problem.  For  the 
dynamic  problem,  then,  more  is  to  be  gained  by  having  m 
large . 

Or,  what  can  be  gained  if  we  allow  two  transfers  per 

period  instead  of  one?  With  a single  transfer  we  have  a 

2 2 

flow-time  (m-^co)  of  b aq(3“1/k)/k  . With  two  transfers 

the  value  is  b^aq{5/2  - 1/2k;/k^.  The  gain  is 

1 - 5/2  - -l/2k  = 1 - 1/k  = 17J 

3 - 1/k  2(3  - 1/k) 

for  k large. 

These  numerical  values  should  not  be  taken  too 
literally.  But  they  do  give  qualitative  insignts  into  the 
relative  optimality  of  suggested  schemes. 


i-Ji  PYOMlfi.  Multlplc-Bua  Problem:  Transfer-Free 

Let  us  lastly  consider  some  dynamic  multiple  bus 
schemes  with  no  transfers.  As  in  the  static  case  of  Section 
2.3.3«  we  do  not  guarantee  optimality.  However,  simple 
Interesting  schemes  can  still  be  studied. 


Two-Bua  Seheaea 
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Conalder  the  following  acheme,  t'(2),  aa  ahown  In 


Figure 

3.1.  Partition 

R 

into  m 

subregions . 

Bus 

I 

visits 

the  subregions  in 

order 

('•1  .'*2. 

...r^j,...). 

and 

bus 

11 

visits  them  in 

order 

( rj, , rij.') 

,...r'],...}. 

Bus 

I 

serves 

demands  of 

the 

form 

(ri,rj) 

where 

I ( J-1 jmodm I < m/2,  and  bus  II  servea  deaanda  with 

l(J-i}aodml  > m/2.  Demands  (ri,rj^}  are  snared  equally 
between  the  buses. 


Figure  3.1  Two -bus  dynomic  scheme  without  transfer 

For  this  Bobeae,  let  6 be  the  tiae  spent  by  each  bus 

in  each  of  its  subregions.  During  a tiae-interval  0,  qO 

demands  (l.e.  2qd  points)  arise.  In  equllibrlua,  each  bus 

Bust  collect  qe  points  in  each  subregion,  l.e. 

9 • by^ya/B , and  9 ■ b qs/B.  The  time  for  a bus  to  make 

2 

a circuit  of  the  b subregions  Is  P > bO  ■ b qa.  For  b 
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large , 

the  avera 

ge  wait-time 

is 

?/2 

and 

the 

average 

travel 

-time  is 

P/4.  Thus, 

for 

m 

large , 

the 

average 

flow-time  for  T* 

(2),  viz. 

„’(2) 
rq  , 

is 

' ( 2 ) 2 

Fg  r 2b  qa  asymptotically  a.e. 


A second  two-bus  scheme  t"(2}  is  similar  to  that 
above,  but  with  both  buses  travelling  through  the  m 
subregions  in  the  same  order  ( r i , r2 , • . • rn , . . . ) . Each  bus 
always  visits  as  many  points  as  possible  in  each  subregion 
it  enters.  We  can  suppose  too  that  the  buses  are  staggered, 
with  bus  II  in  region  ^ ^ )niodm  ^ 


region 

»'i- 

For  this 

scheme 

we  have 

that 

, when 

m 

is 

large , 

the 

average 

travel- 

time  is 

P/2 

and  the 

average 

wait-time  is 

P/4,  where  P 

: b^qa  again. 

Thus  , 

if 

the 

average 

flow 

-time  is 

-"(2) 
r q , 

then  for 

m large, 

" (2)  2 

Fg  s 2b  qa  asymptotically  a.e. 

k 


Although  both  schemes  T*  and  t"  have  the  same 
average  flow-time,  T*  may  be  more  desirable  for  a number 
of  reasons.  First,  for  T*  the  average  wait-time  is  longer 
and  the  average  travel-time  shorter  than  those  of  t".  So, 
with  t',  passengers  have  more  'coffee  time  at  home*. 
Second,  with  t”  the  buses  will  not  always  remain  staggered 
beoause  of  additional  disturbances  and  uncertain 
travel-times.  In  fact,  the  buses  will  tend  to  merge 
together  (a  similar  phenomenon  is  described  by  Newell  [15]) 
so  additional  controls,  causing  additional  delays,  will  be 
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needed.  Third,  suppose  that  R Is  Itself  a 'subregion* 
within  a larger  system,  and  with  a transfer  point  Just 
before  r^,  say.  Then  t"  is  unfair  in  that  a passenger 
originating  in  r^  will  always  have  a longer  travel-time 
than  one  originating  in  rg,.  The  variance  of  the 
travel-times  will  be  smaller  for  T*. 

1 Three-Bus  Scheme 

Finally,  consioer  a three-bus  scheme  TO)  that  is  a 
simple  extension  of  t'(2)  above.  The  region  R is 
aivlded  into  3m  subregions,  each  with  area  a/3m.  The 
buses  perform  routes  as  shown  in  Figure  3.2;  a passenger  is 
served  by  the  most  suitable  bus,  that  is,  one  that  will 
provide  a shortest  travel-time  for  him.  Let  6 be  the  time 
spent  by  a bus  in  each  of  its  subregions.  During  this  time 

2q0  points  arise,  and  for  equilibrium  2qe/3  points  must 

be  visited  by  the  bus.  So,  6 s by2qe/3ya/3m  and 
2 

0 * 2b  qa/9m  asymptotically  a.e.  Each  bus  has  period 
2 

P > 2md  s 4b  qa/9.  When  m is  large,  the  average 

wait-time  for  a passenger  is  P/2.  The  average  travel-time 

ia  clearly  less  than  this,  and  can  be  shown  to  be  P/4. 

Hence  the  average  flow-time  Fq^^  is  3P/4,  l.e.  when  m 

is  large. 


( q ) .2 

Fq  s ” 9*  ssymptotioally  a.e. 


r 
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We  conjecture  that  this  Is  the  minimal  average 
flow-time  for  the  three-bus  transfer-free  case.  However, 
note  again  that  a passenger  requiring  service  from  ri  to 
rgjji  has  a very  long  travel-time  even  though  the  regions 
are  adjacent.  A six-bus  scheme  might  have  two  buses 
travelling  in  opposite  directions  on  each  of  the  routes 
shown  in  Figure  3.2,  and  this  would  reduce  the 
discrimination . 


f CH&ETER  11 

I I.QWARPS  Iii£  PRACTI.CAL  QL  DIAL-A-RIDE  SYSTEMS 


Hi  1 Introduotlon 

In  Chapters  II  and  III  we  have  presented  a mathematical 
theory  for  idealized  bus  problems,  which  we  have  tacitly 
claimed  is  relevant  to  the  real  dial-a-ride  problem.  In 
order  to  Justify  this  claim  we  must  relate  the  real  and 
idealized  problems  and  draw  upon  the  theory  to  obtain 
implications  for  real  design.  In  fact,  we  would  like  to  go 


even  further  and 

propose 

an 

approach  for 

dial-a-ride 

scheduling,  which 

is 

of 

both 

theoretical 

and  practical 

Interest.  It  is  this 

that 

is 

the  goal  of 

the  present 

chapter . 

If  a dial-a-ride  system  is  required  in  a certain  area, 
how  should  the  service  be  designed  and  scheduled?  Our 
response  is  essentially  described  by  the  following  three 
planning  steps.  First,  at  the  most  basic  level,  it  must  be 
decided  what  approach  is  to  be  taken.  The  proposed  approach 
is  through  the  decentralization  of  the  large  problem:  each 
bus  la  required  to  specialize  only  in  passengers  of  a 
certain  'type',  determined  by  their  origin-destination  pair. 
Once  this  approach  has  been  adopted,  the  second  aspect  of 


design  Is  that  of  designating  the  decomposition;  to  do  this 
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we  must  specify  the  roles  to  be  played  by  the  Individual 
buses,  we  must  select  subregions  and  plan  transfers.  At  the 
tnird  level  Is  the  more  aetaiied  operational  design;  here, 
we  must  establish  decision  rules  that  will  determine  the 
behaviour  of  the  buses  when  passengers  arrive. 

This  trinity  yields  a class  of  scheduling  algorithms 
which  we  wish  to  propose.  The  class  is  described 
essentially  by  five  'elements'  in  Section  M . 2 . We  draw  upon 
the  theory  of  Chapters  II  and  III  to  motivate  these 
qualitative  elements  and  also  to  derive  quantitative 
techniques  for  cesign.  In  Section  4.2  too,  some  of  the 
theoretical  assumptions  and  their  limitations  emerge  more 
Clearly.  Adjustments  are  therefore  required  if  the 
quantitative  approach  is  to  be  practically  meaningful,  and 
some  of  these  are  discussed  in  Section  4.3.  Finally,  to 
give  an  Improved  grip  on  the  ideas,  we  describe  an  example 
In  Section  4. 4,  a 'pseudo-real  situation'.  In  which  the 
approach  is  applied. 


»L^  A Approach 

This  section  describes,  in  some  generality,  an  approach 
towards  the  design  of  dial-a-ride  scheduling  algorithms  with 
which  we  Shall  be  concerned.  Given  are  a region  where  the 
transportation  system  Is  to  be  designed,  and  an  anticipated 
probability  distribution  according  to  which  passenger 
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demands  will  arise.  A certain  scheduling  facility  (a 
computer?),  receiving  these  demands,  is  to  determine 
suitable  tours  for  a fleet  of  buses. 

The  following  are  five  elements  that  characterise  a 
class  of  dial-a>ride  scheduling  algorithms: 

(1)  Decomposition  of  Area  and  Specialization  of  Buses 

(11)  Simple  Tours 

(111)  Rules  for  Visiting  Subregions 
(iv)  Transfer  Points 
(v)  Method  of  Quotation. 

Our  interest  In  these  algorithms  Is  motivated  by  the  models 
of  Chapters  II  and  III.  We  shall  discuss  the  five  elements 
In  detail  and  draw  upon  the  theory  to  explain  and  to  Justify 
them.  For  this,  It  must  be  assumed  that  the  Idealized 
models  do  indeed  focus  upon  the  qulntessentlals  of  the  real 
problem.  Of  course,  for  tractablllty , many  assumptions  and 
simplifications  have  had  to  be  made;  some  of  the  more 
important  of  these  will  be  discussed  In  Section  4.3. 
Despite  the  assumptions,  the  many  reasonable  qualitative 
Implications  do  yield  Interesting  algorithms.  And  the 
algorithms,  when  extended  to  the  real  problem,  have  many 
features  of  practical  Interest  to  boot, 

Undoubtedly,  fundamental  to  the  idealized  algorithms  in 
Chapters  IZ  and  III  is  their  decomposed  and  discrete  nature: 
the  buses  specialize  in  certain  demand-types  only,  and  the 
tours  are  all  'simple*.  For  example,  recall  Theorems  2 and 


I 
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3 and  the  tour 

Tn  which 

approximate  to 

any  required 

accuracy  the 

optimal 

performance ^ ^ ^ . 

Even 

very 

simple-minded 

algorithms 

of  this  form 

yield 

good 

performance . 

Thus  we  are 

led  to  claim  that 

the  following 

two  elements  of 

design  are  worthwhile. 

Ci)  Decomposition  Area  ajul  Specializatlfifl 

Buses 


In 

1 tlal 

ly,  th 

e regio 

n in 

which  the  transpo 

rt 

ation 

sy 

stem 

is  to 

be  cons 

true  t 

ed , is  partitioned 

i 

nto  a 

la 

rge 

number 

, ni, 

of 

subregions.  Th 

en 

. a 

pa 

sseng 

er  req 

nesting 

serv 

ice  from  a subregi 

on 

^i 

to 

a subregio 

n rj 

can 

be  directly  cla 

ss 

ified 

ac 

cording  to 

his 

deman 

d-type,  (ri.rj) 

- 

there 

ar 

2 

e m 

such 

demand 

-type 

8.  Now,  assign  t 

0 

each 

of 

the 

buse 

s a fi 

xed  s 

et  of  subregions  w 

hi 

chit 

Vi 

sits . 

Each 

bus  i 

s to 

specialize  within 

its 

su 

bregi 

ons  i 

n serv 

Ing  c 

ertain  demand-type 

s 

only . 

Th 

at  i 

s , it 

visit 

s on 

ly  their  origin 

S, 

or 

de 

stlna 

tions , 

or 

both 

. A bus  may 

also 

Pe 

rlodl 

cally 

visit  • 

trans 

fer  points' . At 

these , 

it 

eolleots 

additional  passengers  and  it  d 

el 

ivers 

a 

subse 

t of  1 

ts  curr 

ent  passengers,  each  o 

f 

whom 

wi 

11  ba 

collaoted  by 

some 

other  bus. 

aeasure 

of 

perforaanoe 

is  distance  trave 

11 

ed  or 

•v«ras«  flow-time.  The  soheaes  have  thus  been  found  to  be 
robust,  alniaizlng  eriterls  relevant  both  to  the  operator  of 
the  aystea  and  to  the  passengers.  We  shall  dlsousa  these, 
and  other  aeaaures  of  perforaanoe,  subsequently. 
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By  the  dl 

scussion 

in 

Sectl 

on  2 . 3 1 

we 

should  assign 

buses 

to  the 

demands 

in 

such  a 

way  that 

all 

buses  share  an 

equal 

portion 

of  the  ' 

load 

' and 

that  each 

bus 

is  busy  all 

the  time.  We  should  also  ensure  that  the  buses  are  always 
suitably  dispersed  over  the  area:  this  will  promise  good 
service  to  the  future  demands  that  arise  according  to  the 
given  probability  distribution. 

The  number  m of  subregions  to  which  each  bus  is 

assigned  is  important,  particularly  in  dynamic  problems. 

However,  this  number  does  not  need  to  be  too  large  before 

good  performance  is  attained  (see  Sections  3.1  and  3>2}. 

When  the  demand  distribution  is  uniform,  the  subregions  can 

( 2 ) 

all  simply  be  of  equal  area,  yet  otherwise  arbitrary  In 

practice,  however,  physical  constraints  such  as  boundaries 
of  towns,  major  roads,  and  so  on,  will  influence  this 
choice . 

Once  the  buses  have  specialized,  the  demand-type  of  a 
passenger  who  requests  service  will  trivially  determine  the 
bus,  or  sequence  of  buses,  to  which  he  is  to  be  assigned. 
The  central  scheduler  that  receives  the  demand  performs  no 
long  search,  but  only  Informs  the  relevant  bus  and  provides 
the  book-keeping. 

HU  SilBlt  IfiiUlA 

Passengers'  demands  arise  over  time.  Depending 


(2)  Additional  remarks  will  be  made  in  Section  4.3 
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upon  the  actual  requeats,  tours  for  the  buses 
within  their  subregions  must  be  deternined.  Any 
point  which  a bus  visits  must  be  'feasible*  - it 
must  be  an  origin  or  a destination  of  a passenger 
who  is  already  aboara.  (Tnese  points  are  of 
course  within  the  bus's  speciality.)  The  points 
are  to  be  visited  in  some  'optimal'  order.  The 
problem  is  simplified  by  constructing  tours  in 
discrete  stages,  whenever  a bus  enters  a new 
subregion.  At  this  time  of  entry,  let  the  bus 
calculate  an  optimal  travelling  salesman  tour  in 
the  upcoming  subregion;  any  point  visited  there  is 
feasible  at  this  Instant  of  time.  The  calculation 
is  to  be  done  on>llne. 

We  have  termed  tours  of  this  form  'simple'.  Recall 
that  the  use  of  these  tours  was  Justified  for  the  static 
problems  of  Chapter  li,  and  it  was  claimed  that  they  were 
suitable  for  the  dynamic  problems  too  (Chapter  III,  page  43, 
footnote  (1)). 

If  the  subregions  are  small  and  only  a few  points  are 
to  be  visited  in  each,  it  may  be  easy  for  a driver  to  'see' 
what  route  he  should  take.  Otherwise,  he  might  have  a 
minicomputer  aboard  that  finds  a tour  for  him.  Efficient 
approximation  algorithms  for  the  travelling  salesman  problem 
are  given  in  [12,3],  In  this  way,  the  oaloulation  of  tours 
is  decentralized  - it  reduces  to  a number  of  small  discrete 
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probleffls.  By  choosing  the  numoer  of  subregions  carefully, 
it  is  possible  to  guarantee  that  the  total  computation  time 
is  bounded  by  a polynomial  function  of  q,  the  arrival-rate 
of  passengers  (see  Karp  [9]).  Note  too  that  the 
communication  between  the  buses  and  the  central  scheduler 
does  not  need  to  be  continuous,  so  relatively  cheap 
communication  equipment  should  suffice. 

(ill)  Rules  lojc.  Yialting  Subreglona 
Two  questions  remain  whenever  a bus  exits  from  a 
subregion  at  a time  ti:  which  region  should  it 
enter  next,  and  which  points  should  it  visit 
there?  Well,  in  any  region  it  enters,  it  must 
visit  all  points  within  its  specialization  that 
were  feasible  there  at  time  And,  it  should 

visit  its  regions  successively ^ ^ \ in  a fixed 
order . 

These  rules  were  proved  valid  for  the  case  in  which 
demands  were  uniformly  distributed,  in  Section  3.1. 
Essentially,  they  require  that  buses  'always  do  as  much  as 
is  possible'.  They  imply  that  a bus  will  not  deviate  in  a 
backward  direction  from  its  prospective  route  to  service 


(3)  This  should  be  further  qualified.  It  holds  when  the 
distribution  of  points  is  uniform.  In  the  ease  of  more 
general  distributions,  a generalized  rule  may  be  required. 
For  example,  visit  next  the  region  with  largest  'point 
density*  (i.e.  the  number  of  feasible  points  per  unit  area) 
- see  Section  4.3(1). 
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newly  arriving  demands,  as  this  may  delay  future  passengers. 
After  t^,  feasible  points  that  arrive  in  the  subregion  are 
not  collected:  a passenger  who  Just  misses  the  bus  must 
wait  for  the  next  one  to  enter  the  subregion.  Even  though 
some  passengers  are  inconvenienced,  tneir  contribution  to 
the  average  performance  is  kept  small. 

It  may  happen  that,  in  a particular  subregion,  there 
are  no  feasible  points.  Then,  the  bus  can  simply  bypass 
tnat  subregion  or  travel  through  it  as  efficiently  as 
possible  to  the  next.  We  can  think  of  the  bus  as  t'ravelllng 
through  Its  subregions  along  a fixed  route  and  deviating 
within  each  subregion  to  collect  and  deliver  passengers 
there.  Passengers  then  know  the  direction  in  which  the  bus 
is  proceeding,  and  should  be  happy,  knowing  roughly  where 
they  are  headed. 

(iv)  Transfer  Points 

The  transfer  points  (see  (i)  above)  at  which  a bus 
stops  from  time-to-tlme  are  typically  on  the 
boundary  of  some  subregions.  They  must  be 
considered  during  the  calculation  of  the  optimal 
travelling  salesman  subtours.  Our  description  has 
been  sufficiently  general  to  allow,  for  example, 
that  some  buses  have  'null  subregions',  travelling 
only  between  transfer  points;  tnese  are  the 
familiar  line-haul  buses. 


r 


-73- 


Transfer  points  were  shown  to  be  necessary  in  Sections 
2.3.2  and  3.3  if  the  delay  incurred  can  be  neglected.  For, 
they  provide  vast  InproveBents  in  performance,  and  this 
improvement  increases  with  the  number  of  transfers. 
However,  too  large  a numoer  of  transfer  points  is  clearly 
impractical,  and  passengers  might  attach  an  additional  cost 
to  the  numoer  of  transfers  they  have  to  make.  Also,  the 
theoretical  approach  of  having  all  buses  meet  at  the  same 
transfer  point  may  not  be  cesirable:-  the  additional  trip 
from  the  region  to  the  transfer  point  will  not  be  negligible 
when  the  numoer  of  passengers  is  small.  In  this  case,  the 
line-haul  bus  effectively  creates  the  single  transfer  in  a 
practical  way. 

1y1  Quotation  Times 

In  the  statement  of  the  dial-a-rlde  problem  in 
Section  1.1  it  was  required  that,  at  the  time  a 
passenger  requests  service,  he  be  quoted  a time  of 
collection  and  delivery.  The  scheduling  procedure 
we  have  described  has  not  needed  such  a quotation. 
We  regard  it  as  a separate  secondary  problem, 
requiring  that  an  estimate  of  the  'load*  on  the 
system,  and  hence  an  estimate  of  the  speed  of  the 
buses,  be  made. 
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The  approach  which  we  have  described  guarantees  a 
certain  average  service  to  the  set  of  all  passengers.  Thus, 
if  quotations  are  made  to  individual  passengers,  and  if 
these  are  then  Introduced  as  additional  constraints,  the 
average  service  to  passengers  will  be  reduced.  In  practice, 
quotations  are  required,  but  the  scheduler  need  not  be 
committed  to  them.  We  believe  that  simple  constant 
estimates  of  the  times  of  collection  and  delivery  will 
usually  suffice.  It  is  known  through  which  subregions  a bus 
is  to  travel:  the  time  spent  deviating  in  each  subregion 
could  be  estimated  and  quotation  times  could  easily  be 
provided.  So,  an  accurate  prediction  of  travel-time  between 
points  is  avoided. 

The  elements  (l)-(v}  above  outline,  in  somewhat  general 
terms,  a class  of  scheduling  algorithms.  The  Dasic  approach 
is  described  by  (1)  and  (11)  and  we  have  claimed  that  this 
approach  is  Justified  In  an  asymptotic  probabilistic  sense. 
The  detailed  operational  design  is  described  by  (ill)  for  a 
special  case,  and  in  the  next  section  It  will  be  seen  that 
this  can  be  extended  fairly  generally.  What  is  lacking, 
then,  is  a methodology  for  selecting  a suitable 
decomposition  and  for  organizing  the  manner  in  which  the 
buses  specialize. 

For  the  special-case  idealized  models  of  Chapters  II 
and  III  no  such  methodology  is  required  since  the  solution 
is  relatively  easy:  when  the  demands  are  uniformly 


1 

I 

I 

I 
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diBtrlbuted  over  the  region,  partition  it  into  aubregiona  of 
equal  area.  Then,  in  the  manner  of  Algorithm  6.2  the  buaea 
apecialize  in  diajoint  aubregiona  and  have  cyclic  toura 
through  them  when  the  arrival-rate  of  paaaengera  la  large. 
In  general,  however,  the  anawer  la  not  aa  eaay  and 
dif ficultiea , beaidea  the  nonuniformity  of  dlatrlbutiona , 
abound  in  practice.  Indeed,  let  ua  remark  upon  a few 
practical  conalderationa  which  Indicate  that  a general 
methodology  for  the  determination  of  aubregiona  and 
bua-routea  would  aurely  be  very  complex. 

Flrat,  what  criteria  actually  determine  good  aervice? 
The  criteria  aaaumed  in  Chaptera  II  and  III,  viz.  the 
dlatance  travelled  by  the  buaea  and  the  average  flow-time  of 
paaaengera  are  not  all-important,  but  are  only  flrat 
approxlmatlona  to  real  requi rementa . Tnere  are,  for 
example,  the  variance  of  paaaengera*  travel-tlmea 
(meaauring,  in  a aenae,  the  'fairneaa*  to  the  paaaengera) 
and  the  maximum  of  all  paaaengera*  travel-tlmea  (aince  we 
wlah  to  provide  all  demand-typea  with  good  aervice).  Or, 
there  are  queationa  of  the  reliability  of  the  ayatem  (what 
happena  if  a bua  breaka  down?)  and  ita  atablllty  (what 
effect  will  additional  delaya  and  uncertain  travel-time 
have?) . 

Furthermore,  a new  dial-a-rlde  ayatem  would  have  to 
suit  the  existing  oiroumstanoas.  Boundaries  of  towns, 
rivers  or  roads  may  impose  physloal  oonstralnts  on  the  way 
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in  which  the  region  can  be  subdivided.  We  may  wish  to 
supply  a special  quality  of  service  to  certain  demand-types 
le.g.  to  disabled  persons,  or  see  Example  <1.4).  A 
dial-a-ride  system  must  be  compatable  with  existing  public 
transportation  systems:  it  should  supplement  them,  and  yet 
could  also  depend  upon  them.  For  example,  existing  bus-  or 
subway-routes  might  reduce  the  number  of  line-haul  buses 
that  are  needed  if  the  subregions  are  designed  to  take 
advantage  of  them. 

Thus,  dlal-a-rlde  schemes  must  be  tailor-made  to  suit 
inoividual  complex  requirements.  Even  though  a general 
methodology  has  not  been  proposed,  our  study  has  yielded 
many  qualitative  insights,  and  no  doubt  additional 
principles  could  be  obtained  with  further  research. 

On  the  other  more  favourable  hand,  however,  we  wish  to 
suggest  that  our  approach  does  yield  a practical  design 
tool.  The  emphasis  throughout  Chapters  II  and  III  was  upon 
the  comparison  of  algorithmic  performance  (see  Definition 
(2.1)}.  In  general,  it  is  possible  to  quantify  the 
performance  of  a suggested  scheme  whenever  a distribution  of 
demands,  a decomposition  of  the  region  into  subregions  and 
an  organization  of  specialist  buses,  are  given.  For,  one 
can  estimate  the  time  spent  by  each  bus  in  each  of  its 
subregions  (using  Beardwood's  formula.  Theorem  1)  and 
evaluate  its  dynamic  behaviour  as  it  travels  around  its 


-77- 


( 4 } 

subregions  and  visits  Its  transfer  points  . The  service 
supplied  to  the  various  deoand-types  by  the  suggested  scheme 
can  tnen  be  evaluated  with  respect  to  average  flow-time  and 
to  many  other  performance  criteria.  Thus,  a suggested 
scheme  can  be  quantitatively  studied.  Usually,  closed-form 
analytic  formulas  should  be  obtainable  to  express  Its 
cnaracteristlcs . Then,  Investigations  can  be  made  to 
analyse  the  behaviour  of  the  scheme  when  the  parameters 
e.g.  the  number  of  buses,  the  demand  distribution,  the 
number  of  transfer  points,  and  so  on  - are  varied.  Finally, 
alternative  schemes,  with  alternative  regional  decomposition 
and  bus-speclallzatlon  can  be  evaluated  and  compared.  These 
are  analytic  studies,  and  avoid  the  need  for  complex 
simulations.  An  example  that  briefly  Illustrates  such  a 
quantitative  investigation  is  given  In  Section  4.4. 

It  Is  In  order,  finally,  to  make  some  additional 
remarks  on  the  design  approach  which  we  have  expounded  In 
this  section. 

The  approach  can  be  applied  on  a broader  scale.  The 
'city*  may  Itself  be  a module,  connected  to  others  within  a 
larger  system.  Also,  the  method  Is  relevant  to 

'one-to-many*  systems,  [21],  where  passengers  have  a common 
origin  and  wish  to  be  taken  to  various  destinations.  For, 
It  is  besy  to  assign  Inoomlng  passengers  to  buses 


(4)  This  was  illustrated  for  the  simple  models  in  Chapter 
XXX. 
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simultaneously , merely  treating  the  single  origin  as  a 
transfer  point.  ' Many-to-one * systems  are  similar. 

In  [21],  the  existence  of  'advance  requests',  where 
passengers  request  service  some  time  in  advance,  is  regarded 
as  important.  With  our  approach  these  demands  are  easy  to 
handle.  If  the  'period'  of  a bus  is  P - this  is  the  time 
it  takes  for  the  bus  to  make  a circuit  of  all  its  subregions 
the  passenger  can  be  promised  collection  during 
[ti.t^-fP]  for  any  suitable  time  t^.  Then,  when  the 

region  is  entered,  the  passenger  will  be  collected  if  this 
can  be  during  [t^.t^-fP],  otherwise  the  collection  is 
deferred.  Delivery  by  a certain  time  can  similarly  be 

guaranteed  in  advance. 

There  is  a further  important  implication  of  the 
decomposition  described  in  (1)  above.  The  success  of  a 
transportation  system  depends  ultimately  upon  the 
co-operation  of  the  bus-drivers,  particularly  if  they  have 
decision-making  capabilities.  Since  each  bus  has  certain 
fixed  demand-types  assigned  to  it,  it  is  possible  to  measure 
the  efficiency  of  individual  drivers.  If  each  driver 
travels  as  well  as  possible  through  his  subregions,  the 
total  system  performance  will  be  maximized.  So,  some 
inoentlve  might  be  given  to  the  drivers  that  will  encourage 
them  to  perform  well  and,  for  example,  to  reduce  the  number 
of  unnecessary  stops.  Furthermore,  since  the  drivers 
communicate  with  the  central  scheduler  only  intermittently. 


I 
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the  ' big-orother ’ paranoia  of  continuously-watched  drivers 
is  reduced. 


iL^  OuantitatlYfi  Aagecta 

In  Section  4.2  it  was  described  how  suggested  schemes 
can  be  investigated  by  quantifying  their  behaviour. 
Investigations  such  as  tnese  are  limited  by  the  many 
theoretical  assumptions  of  Chapters  II  and  III,  since 
inaccuracies  are  then  produced  in  the  resulting 
measurements.  In  this  section  we  discuss  some  of  the  more 
important  of  tnese  assumptions  and  investigate  to  what 
extent  they  can  be  relaxed.  We  Indicate  how  the 
inaccuracies  might  be  reduced  by  making  adjustments  in  the 
applied  formulas. 

1 The  Demand  Distribution 

In  deriving  formulas  measuring,  for  example,  the 
average  flow-time  for  a suggested  scheme  (see  Theorems  5 and 
6)  Tneorem  1 is  used  to  obtain  the  time  spent  by  a bus  in 
each  of  its  subregions.  For  this,  it  must  be  assumed  that 
the  points  which  are  visited  In  each  subregion  are  drawn 
from  a uniform  probability  distribution  there:  this  was  the 
case  in  our  previous  analyses  since  the  demand  pairs  were 


' random' 


W ti  A V k4 


transportation  systen  Is  to  be  designed,  and  an  anticipated 
probability  distribution  according  to  which  passenger 
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Suppose,  now,  tnat  the  points  to  be  visited  in  a 
subregion  are  not  uniformly  distributed,  but  arise  according 
to  a probanility  density  function,  p,  say.  Then  it  has 
been  shown  in  [2]  that  the  optimal  tour-length  Ln 

satisfies 

lim  Ln  = b/  p^^^dv  a.e. 

n->oo  Vn  R 

When  p s 1/a,  we  have  the  uniform  distribution.  Since 
1/2 

f p dv  <,  yS,  nonuniformity  of  the  distribution  decreases 

R 

Ln*  Consequently,  the  additional  structure  in  the 

nonuniform  distribution  can  be  used  to  advantage. 

The  expression  above  might  be  used  to  improve  the 
analyses  when  the  distribution  is  nonuniform,  and  when  all 
origins  and  destinations  are  inoependent  of  one  another. 
The  discretization  imposed  by  simple  tours  will  remain 

• useful,  however.  For,  if  the  whole  region  is  diviced  into 

j 

m subregions  and  if  we  approximate  the  distribution  p to 
be  uniform  within  each  subregion,  then  by  making  m large, 
the  approximation  can  be  made  arbitrarily  negligible. 

» 

Now,  consider  a bus  visiting  its  subregions.  What  rule 
(l.e.  element  (ill])  should  it  use?  It  is  not  hard  to  see 
test  in  order  to  alnlmlze  average  flow-time  the  bus  should 
always  *do  as  auoh  as  possible’  (compare  with  Lemma  5.1). 
Taaa.  waeaaver  it  enters  a subregion  it  should  always  visit 
I ••  saay  paints  as  possible  there.  But,  which  subregion 

.»«.  • It  visit?  If,  In  each  of  its  subregions 


o^  fsaalble  points  to  be  visited 
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and  If  then  Ti(nj^)  will  be  the  time  spent  In  each,  choose 
1 so  as  to  mlnlDlze  T ( n^^ ) /n  . Even  when  the 
arrival-rate  q of  demands  Is  large,  It  may  be  tne  case 
that  the  bus  does  not  visit  Its  subregions  successively.  In 
general,  the  analytic  calculation  of  average  flow-times  will 
be  more  difficult  than  before,  but  will  still  be  feasible. 

Note  that  these  remarks  are  not  true  if  the  origin  and 
destination  in  a demand  pair  are  correlated. 

Finally,  It  Is  in  order  to  comment  upon  the  stochastic 
process  by  which  demands  arrive  over  time.  We  assumed  a 
constant  arrival-rate  q,  and  the  analysis  could  be 
extended  for  a Poisson  process  with  mean  q.  In  practice, 
q will  change  over  time.  This  change  will  be  relatively 
slow  - at  certain  times  of  the  day  alternative  designs  for 
the  decomposition  and  specialization  might  be  selected,  so 
as  to  suit  the  changing  operating  conditions.  It  Is 
believed  that  a few  discrete  such  changes  should  suffice  in 
practice . 

-0-^.2  ^ Metric 

We  have  taken  as  the  distance  between  any  two  points  In 
the  city  their  euclidean  separation  - i.e.  the  length  of 
the  straight  line  Joining  them.  In  reality,  the  distance 
travelled  by  a bus  is  not  as  simple.  If  the  actual  distance 
can  be  approximated  by  a constant  multiple  of  the  euclidean 
distance,  the  analysis  of  Chapters  II  and  III  remains  valid 


(1)  Th«  ■•■sure  of  perforaanoe  Is  distance  travelled  or 
average  flow-tlae.  Tbe  aobeBes  have  thus  been  found  to  be 
robust,  ■inlaiaing  criteria  relevant  both  to  tbe  operator  of 
tbe  ayatea  and  to  the  passengers.  We  shall  discuss  these, 
and  other  aesaurea  of  perforaanoe,  subsequently. 


a 
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wlth  very  minor  modifications. 

Alternatively,  the  distance  between  two  points  x,  y 
may  be  better  approximated  by  a multiple  of  the  'Manhattan 
metric',  or  l\  given  by 

distai.ce(x,y)  = Ixi-Xll  + \^2~y2^  • 

Again  the  results  hold,  this  time  with  a different  value  of 
the  constant  b [2]. 

More  general  metrics  that  take  into  consideration 
Indirect  or  faster  roads,  promise  to  be  an  order  of 
magnitude  more  difficult;  for,  there  is  then  no  result 
corresponding  to  Theorem  1. 

3t.3  Aayaptgtlc  Approximation 

Possibly  the  most  glaring  Inadequacy  of  the  approach  is 
that  it  is  asymptotic  in  the  number  of  passengers.  The 
expressions  which  are  obtained  are  laws  of  large  numbers, 
and  require  that  the  number  of  passengers  and  their 
arrival-rate  q be  large.  To  apply  the  formulas  we  would 
Ideally  seek  a 'central  limit'  result  allowing  us  to 
estlaate  the  probable  magnitude  of  tbe  errors  when  q is 
saall,  and  to  deteralne  the  size  of  q needed  for  reliable 
estiaates . 

Sadly,  a general  result  of  this  nature  is  not 
available.  A theoretical  analysis  is  hard  and  results 
depend,  for  ezaaple,  upon  the  shape  of  the  region.  For  R 
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square  , 


analytically 

that^^^ 


the 

on 


oest  that  It 
the  convergence 


has  oeen 
of  Theorem 


possiDle  to  do 
1,  Is  to  show 


^ ^ ^n  ^ + OClogn/n) 

Vn 

where  ECL„]  Is  the  expected  length  of  an  optimal  travelling 
salesman  tour  for  a random  Instance  of  the  problem  of  size 
n.  It  might  be  hoped  that  the  actual  rate  of  convergence 
is  faster  than  this. 


In  this  spirit,  let  us  investigate  empirically  the 
convergence  of  Theorem  1.  We  shall  consider  only  points 
which  are  distributed  randomly  in  a unit  square,  and  shall 
seek  to  estimate  E[Lq]  as  a function  of  n.  Comparing 
E[L|]]  with  bVn  will  yield  an  expected  error  In  the 
prediction.  This  error  will  enable  us  to  bound,  in  a sense, 
the  error  due  to  the  asymptotic  approximation  In  dynamic 
versions  of  the  problem  as  well. 

The  value  of  was  estimated  for  ns20,30,...  as 

follows.  Random  travelling  salesman  problems  were  generated 
by  scattering  n points  on  the  unit  square  (x  and  y 
coordinates  were  generated  uniformly  and  Independently).  A 
number  of  Instances  of  each  problem  size  n were  solved^^^. 

(5)  This  can  be  derived  in  the  manner  of  Beardwood'a  Lemma 

6. 

(6)  It  should  be  remarked  that  a similar  experiment  was 
performed  by  Chrlstofldes  and  Ellon  [4j;  however,  details  of 
their  numerical  results  and  the  algorithms  used  were  not 
provided . 


ar«  given  In  [12,31.  In  this  way,  tbs  calculation  of  tours 
Is  dsosntrallzed  - It  reduces  to  a nuaber  of  snail  discrete 
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The  algorithm  used  was  a heuristic  of  Lin  and  Kernighan 
[12]:  it  uses  modest  computation  time  and  has  been  found  to 

'practically  guarantee  optimality'.  For  n large  the 
computation  time  was  still  high,  however,  and  fewer  problems 
were  sampled.  The  results  obtained  are  tabulated  in  Figure 
4.1.  From  these,  we  estimate  the  value  of  b to  be  about 
.765;  Figure  4.2  compares  the  empirical  expected 

tour-length  with  the  value  of  .76575.  In  Figure  4.3  we 

plot  E[Ln]/7n  versus  n,  indicating  how  it  converges  to 
.765.  The  curve  .765  + 4/n  is  given  for  comparison:  it 

seems  clear  that  E[L„]/7n  converges  at  a faster  rate  than 
0(1/n). 

It  is  interesting  to  note  that  for  every  random 
problem-instance  tried,  with  b * .765, 
b7n  i Lfl. 

So,  with  high  probability  the  tour-length  predicted  by 
Tneorem  1 is  leas  than  the  actual  value.  Indeed,  from 
Figure  4.3, 

bTn  i E[Ln]  i bTn  ♦ 4/n  (4.1) 

This  yields  a fairly  close  bound  on  the  magnitude  of  the 
error  Implicit  In  Theorem  1. 

Now  let  us  obtain  an  estimate  of  the  convergence  rate 
for  a bus  problem  - we  arc  interested  primarily  In  the 
dynamic  version.  In  evaluating  the  performance  of  a scheme 
the  time  spent  by  a bus  In  each  of  Its  subregions  Is 
estimated,  so  we  can  utilize  the  results  above.  We  consider 


i 


CHAV 

For  txaapl*,  viait  next  the  region  with  largest  'point 
density*  (i.e.  the  nuaber  of  feasible  points  per  unit  area)  f 

- see  Section  4.3(1). 
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Asymptotic  convergence  to  Beordwood's  tormulo.  For  each  n,  the  range  of  experimental 
results  is  shownj  ECLpI  is  denoted  by  x.  For  comporison,  the  curve  .765/r”  is  plotted. 
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only  a very  simple  problem:  a single  bus  Is  serving  a 
region  of  unit  area  and  visits  m subregions  successively. 
Each  subregion  la  a square  with  area  1/m,  and  the  system 
Is  In  steady-state. 


Let  q be  the  arrival-rate  of  demands,  and  let  E[Pq] 

be  the  expected  period  - l.e.  the  time  for  the  bus  to  visit 

all  m subregions.  In  Equation  (3.4)  we  predicted  that 

2 

ElPq]  would  behave  asymptotically  as  2b  q.  We  now  ask 
the  following  simplified  question:  can  we  bound  the 

relative  error  In  this  prediction  by  a function  of  q.  That 
la,  seek  a bound  to 

E[Pn]  - 2b^q 
2b^q 

Let  us  make  another  simplifying  assumption.  Suppose 

that,  as  the  bus  travels  through  Its  subregions.  It  visits 

(7 ) 

the  same  number  of  points,  n.  In  each  . Then,  the 

expected  time  spent  In  each  subregion  Is  E[Ln3'/T7m,  and 
the  expected  period  Is  E[Pq]  x mE[Lu]y1 /m.  During  a 
period,  mn  points  are  visited, 

l.e.  q X n (^*.3) 

2E[Ln3'/^ 

(7)  This  assumption  drastically  simplifies  the  analysis. 
The  number  of  points  to  be  visited  In  each  subregion  Is 
really  random  and  we  get 

00 

Expected  time  In  a subregion  > H pq(n)E[LQ]'/T7m 

nxO 

where  PqCn)  Is  the  probability  that  there  are  n points 
In  the  subregion  when  the  arrival-rate  Is  q.  (Note  that 
there  Is  a correlation  between  the  number  of  points  n 
visited  In  successive  subregions:  there  will  be  sequences 
of  large  n and  sequences  of  low  n.) 


(4.2) 
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From  (^.2), 


error 


Using  (4.1), 


E[Ln3  - b^n/E[Ln] 
b^n/ELLn] 


error  i.  1 (byE  4/n)^  - b^n  I 
' b^n  ' 

s I 8byn  ♦ Ifa/n  | 

w2  2 

b n 

With  b * .765,  the  relative  error  given  by  (4.2)  is  less 
than  12f  when  n,  the  number  of  points  visited  in  each 
subregion,  is  20;  it  is  less  than  when  n is  30; 

and  it  is  less  than  5%  when  n Is  35.  For  any  n, 
(4.3)  determines  the  corresponding  arrival-rate,  q. 


4.1.4  Additional  Uncdrtala.tlaa 

The  following  are  a few  additional  Issues  (mainly 
related  to  stochastic  variations)  that  we  have  neglected  and 
that  we  now  wish  to  note. 

(i)  The  speed  of  a bua  is  not  constant  as  has  been 
assumed.  The  travel-time  between  two  points  is  uncertain 
and  Is  subject  to  traffic  and  road  conditions. 

(li)  Delays  that  occur  when  a stop  is  made  have  been 
negleoted.  A bus  may  have  to  wait  for  a passenger  to  appear 
and  for  his  impedimenta  to  be  loaded.  Similar  delays  occur 
for  alighting  paaaengera. 


(iii)  A bua  is  also  detained  at  transfer  points.  The  tine 
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taken  Is  not  Instantaneous,  but  Is  a function  of  the  number 
of  passengers  who  transfer.  Indeed,  if  buses  are  to  meet  at 
transfer  points  as  was  required  in  Chapters  II  and  III, 
waits  are  Inevitable. 

llv)  From  a passenger's  viewpoint,  the  time  spent  on  a 
lins'^haul  bus  (Including  his  wait}  may  be  appreciable.  In 
our  asymptotic  analysis  these,  as  well  as  the  time  taken  for 
a bus  to  travel  a fixed  bounded  distance  were  neglected. 

(v)  Consider  a bus  travelling  a circuit  around  its 
subregions.  In  Chapter  III  we  calculated  the  equilibrium 
value  of  the  average  flow-time.  If,  now,  the  bus  la 
delayed,  additional  passengers  still  arrive,  and  the  time 
spent  in  the  following  subregions  will  be  longer  than 
anticipated.  So,  the  delay  will  persist  and  the  equilibrium 
will  only  be  re-aehleved  gradually. 

Because  of  these  assumptions,  there  will  be  yet  further 
errors  in  the  derived  formulas.  However,  it  is  believed 
that  many  of  them  could  be  taken  into  account  in  a practical 
analysis.  For  example,  studies  of  urban  traffic  behaviour 
might  be  used  to  more  fully  understand  the  travel-times;  and 
we  might  be  able  to  estimate  the  time  spent  by  passengers  on 
the  line-haul  system.  Then,  the  tools  which  we  have 
proposed  oould  be  further  sharpened  and  calibrated. 
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Finally,  a more  sopnlstloated  technique  for  the 
detailed  operational  control  > one  that  will  take  into 
consideration  the  random  fluctuations  - can  he  suggested. 

Consider  specifying  in  advance  a schedule  of  times  for  the 
buses  to  be  at  transfer  points.  Then,  if  a bus  is  running 
late,  it  can  improve  its  speed  by  selecting  carefully  which 
of  its  passengers  it  should  collect:  some  passengers  remain 
uncollected  until  the  bus  (or  the  next  suitable  bus) 
re-enters  their  subregion.  This  feature  would  require 
additional  machinery  if  it  is  to  be  analysed.  A measure  of 
the  'state'  of  each  bus  would  have  to  be  cefined,  and  it 
would  be  nontrivial  to  determine  the  schedule  in  advance  - 
certain  'slack'  would  be  required  to  allow  for  random 
disturbances.  The  schedule  would  have  to  be  'stable'  with 
respect  to  reasonable  perturbations  and  'flexible',  so  that 
changes  can  easily  be  made  if  things  go  awry.  To  study 
these  topics  would  require  new  methods  of  some  'stochastic 
scheduling  theory'. 

iLJL  An  Example 

Let  us  suppose  that  we  have  a city,  R,  in  which  we 
want  to  construct  a dial-a-ride  transportation  system.  This 
city  is  essentially  composed  of  two  distinct  regions  as 
shown  in  Figure  M.4  - downtown,  R(j,  with  area  a^,  and 
suburbs,  Rg,  with  total  area  ag.  There  is  already  an  j 

efficient  public  transportation  system  operating  In  R^.  ^ 


or  unnecessary  stops.  Purthermore,  since  the  drivers 
ooaaunlcste  with  the  oentrsl  scheduler  only  Internlttently , 


Figure  4.4  An  example  : the  city 

For  most  of  the  day  there  is  heavy  traffic  travelling 
between  Rg  and  Rd  in  both  directions.  The  dial-a-ride 
scheme  is  to  be  designed  for  this  scenario. 

Suppose  we  can  approximate  the  distribution  of  demands 
that  are  to  be  served  by  our  system  by  the  following.  Any 
particular  demand  pair  (o,d}  is  of  the  form  (RgiRd)  with 
probability  p/2,  p6(0,i;;  it  is  of  the  form  (RdtRs^  with 
probability  p/2;  and  it  is  of  the  form  (Rs>^s)  with 
probability  (1-p).  This  distribution,  then,  is  nonuniform. 
The  system  is  not  required  to  serve  demands  (Rd,Rd}' 
Further,  the  value  of  p is  large,  l.e.  fairly  close  to 
unity. 

As  remarked,  the  approach  which  we  wish  to  adopt  is  one 
in  which  R is  decomposed  into  subregions  and  in  which 
buses  spaoiallse.  The  theory  of  Chapters  11  and  111  has  not 
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given  a general  metnoaological  procedure  for  designing  this 
decomposition  when  demands  are  nonuniform,  but  the  following 
could  well  be  reasonable  for  the  present  example. 

In  order  to  discourage  persons  who  travel  between  Rg 
and  R(]  from  using  tneir  private  venicles,  let  us  provide 
them  with  attractive  transportation,  say  a transfer-free 
service.  We  describe  first  a scheme  for  only  these 
(Rai^d)  (I^d'^s^  demands.  Consider  a division  of 

Rg  into  k regions,  r^,...r^,  each  of  area  a^/k,  for 
an  integer  k.  It  might  be  desirable  to  have  natural 
divisions  here,  described  by  boundaries  of  suburbs,  roads 
and  so  on.  This  is  depicted  in  Figure  4.5. 


Figure  4.5  Portition  of  suburbs,  Rj 

Now  allocate  a bus  to  the  region  r^,  for 

ls1,...k;  bus  bj[  will  serve  passengers  with  origin  or 
destination  in  rj,  and  will  travel  also  in  R,).  Each  bus 
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will  travel  a 'simple  tour'  (see  Section  2.2,2),  so 
subdivide  r^,  1=1, ...k  and  R^j  further,  into  m 

subregions  each.  Bus  bj^  travels  a route  as  shown  in 
Figure  4.6.  It  makes  a circuit  around  its  2m  subregions, 
visiting  them  successively.  If  it  enters  a subregion 
at  a time  t,  it  visits  there  all  possible  points  that  are 
feasible  at  t.  These  are  the  origins  of  demands  that  have 
arisen  during  the  time  since  >^i,J  was  last  entered,  and 
all  destinations  in  *'i,J  passengers  who  are  already 

aboard  b^.  The  tour  in  >'i,J  optimal  travelling 

salesman  tour,  and  can  be  calculated  at  t.  If  the  number 
of  points  to  be  visited  in  ^iij  small,  the  driver  might 

easily  be  able  to  choose  his  best  tour.  Otherwise,  there 
might  be  a minicomputer  aboard  that  calculates  the  tour  for 
him.  Efficient  approximation  algorithms  that  might  be  used 
for  this  are  those  in  [12,3]. 


Figure  4.6  Route  token  by  o single  bus  in  o region 
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The  buses  communicate  with  a central  computer  (or 
scheduler)  that  Informs  them  of  the  points  to  be  visited  in 
tneir  subregions.  This  computer  monitors  all  demands  and 
assigns  them  to  the  buses;  the  assignment  decisions  are 
trivial  - no  long  or  hard  search  Is  needed. 


With  k regions  of  Rg  and  k buses,  each  travelling 
at  unit  speed,  the  asymptotic  equilibrium  average  flow-time 
for  the  (Rd,Rs)  and  (Rs,R(j)  demands  in  the  above  scheme, 
can  be  calculated  using  the  techniques  developed  in  Chapter 
III  as  follows.  Consider  the  route  of  a particular  bus. 
Let  0 be  Its  period,  the  time  for  It  to  make  a circuit  of 
its  2m  subregions^  . Let  6d  be  the  time  spent  in  each 
subregion  of  Rdi  and  let  Og  be  the  time  spent  in  each 
subregion  of  rj^.  Then, 

0 s mOg  •*-  mOd  (4.5) 

During  0,  there  arrive  2pqe/k  points  that  are  to  be 
served  by  this  bus,  i.e.  pqO/km  in  each  subregion.  Thus, 
by  Theorem  1 , 


Prom  (*1.5), 


l*».6) 

(*♦.7) 


As  In  the  proof  of  Theorem  5 we  get  tnat  the  average 


(8)  Note  that  we  are  being  rather  loose  here.  These  are  the 
equilibrium  values  of  0 when  q is  large;  the  equations 
to  follow  bold  asymptotically,  with  probability  1. 
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flow-time  of  tnese  passengers  is 

Fn  t)  pqryi:;  ♦ [in^(U1/2m)  (4.8) 

k L J kj 

for  q large.  This  approximation  to  the  value  of  Fq 
should  not  be  taken  too  literally:  some  of  its  limitations 
are  discussed  in  Section  4.3. 

Various  other  similar  schemes  might  also  be  suggested 
for  serving  the  (RsiR^)  and  (RqiRs)  demands.  Suppose 
tnat  the  number  of  buses  is  increased.  Then,  the  average 
flow-time  will  be  minimized  if  we  increase  also  the  number 

of  regions  of  Rg,  and  retain  a single  bus  in  each. 

\ 

However,  it  might  also  be  desirable  to  have  more  than  one 
bus  in  each  region  rj,.  For  example,  two  buses  travelling 
through  the  subregions  as  shown  in  Figure  4.7  would  provide 
a more  equitable  service  to  passengers  (i.e.  » 'r 

variance  of  travel-time)  and  a lower  value  of  the  ma  sum 

flow-time.  This  scheme  is  also  more  reliable:  if  a bus 

breaks  down,  its  passengers  can  still  be  collected  by  the 
otner  bus.  These  are  all  considerations  that  are  important 
in  practice. 

Having  described  how  passengers  travelling  between  R(] 
and  Rg  might  be  served,  we  now  turn  to  suburban  demands, 
of  the  form  (RgiRg).  In  general  we  have  seen  that  a 

sarvioa  with  transfers  will  yield  an  average  flow-time  that 
is  batter  than  that  with  a transfer-free  service.  So, 

consider  a scheme  with  transfers. 
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Figure  4.7  Route  token  by  a pair  of  buses  in  a 
region 

With  J buses,  partition  R5  into  J regions 
r^,...rj  as  was  done  before,  and  assign  a bus  to  each. 
Also  specify  a transfer  point  in  each  region  and  a route  for 
a line>haul  bus,  travelling  between  these  transfer  points, 
in  the  Banner  of  the  'Michigan  Scheae',  as  described  in 
Section  1.1.  This  line-haul  bus  will  play  the  role  of  a 
transfer  point:  the  distance  which  it  travels  is 

asyaptotioally  negligible,  and  it  is  clearly  not  always 
practical  to  have  all  regional  buses  aeetlng  at  the  saae 
transfer  point  as  was  required  in  Algoritha  3<  For  our 
exaaple,  this  is  sketched  in  Figure  4.8. 

Once  aore,  subdivide  the  regions  r ; the  regional 
buses  perfora  siaple  tours  in  their  subregions  in  the  saae 
aanner  as  those  before,  visiting  the  transfer  point  each 
period.  The  previous  reaarks  on  the  nuaber  of  buses  in  eaoh 
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Figure  4.8  Decomposition  for  suburban  demands 


region  are  again  pertinent. 


If  there  are  J regional  buses  and  J regions,  with 
m subregions  each,  the  asymptotic  average  flow-time  can  be 


calculated  from  Theorem  6 (set  ksj;  hsi)  to  be 
2 

Fq  b (l-p)qaa  fa  ♦ j.  - l(1-1/m)*)  ♦ Qi 

? L » J J 


(‘•.9) 


Here,  Is  the  average  time  spent  by  a passenger  on  a 
fixed-route  bus  - It  depends  for  example  upon  the  number  of 
these  buses  and  upon  the  actual  route  taken. 


This  scheme  corresponds  to  a single  transfer  per 
period.  It  Is  possible  to  obtain  further  Improvement  In 
average  flow-time  by  having,  say,  two  transfers  per  period. 
Than,  we  need  2J  regions  of  R^,  with  each  bus  assigned 
to  two  of  them,  and  a transfer  point  In  each.  A bus  visits 
its  two  regions  sltsrnstingly , performing  a tour  through  the 
subregions  in  each,  and  visiting  the  transfer  point 
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tnereafter  (In  the  manner  of  the  tour  of  Algorithm 

C 9 ID 

6).  The  asymptotic  average  flow-time  can  be  calculated  as 

before  to  be 
2 

Fq  b (1-p)qa3  r2.5  ♦ j_  - _L(l-1/m)1  ♦ Q2  (4.9a) 

? L 2m  2J  J 

With  this  approach,  then,  it  is  possible  to  calculate 
the  asymptotic  flow-time  and  other  performance  criteria  for 
various  suggested  schemes  as  in  (4.8)  and  (4.9).  It  is 
indicated  in  Section  4.3  how  adjustments  in  these  formulas 
might  be  made  to  reduce  their  error.  Then,  a parametric 
study  can  be  made  to  analyse  the  behaviour  of  alternative 
system  designs.  For  example,  given  a fixed  number  of  buses, 
determine  how  they  should  be  allocated  between  the  suburban 
and  central  subsystems.  Or,  calculate  the  possible  gain  in 
flow-time  if  the  number  of  buses  is  increased  or  if  there 
are  more  transfers  on  the  suburban  subsystem.  Different 
schemes  with  alternative  decomposition  of  the  area  might 
also  be  studied,  and  other  performance  criteria  such  as  the 
average  waiting-time  or  maximum  travel-time  might  be 
investigated.  This  can  all  be  done  analytically,  avoiding 
the  need  for  complex  simulations. 


CHAPTER  1 


CQNCLUSIQK 


Our  first  aim  In  this  report  has  been  towards  the 
development  of  dial-a-rlde  scheduling  algorithms.  We  have 
attempted  an  analysis  of  the  problem  at  its  most  fundamental 
level,  i.e.  to  question  what  design  approach  should  be 
taken.  In  our  study  we  have  viewed  the  system 
macroscoplcally  - by  considering  the  statistics  of  the 
arriving  demands,  we  have  analysed  probabilistically  the 
combinatorial  aspects  of  the  problem  for  suitable 
mathematical  models.  This  has  led  to  a derivation  of  a 
class  of  algorithms  (or,  schemes)  for  which  performance  can 
easily  be  evaluated,  and  hence  to  an  approach  towards  the 
design  of  dial-a-ride  systems. 

The  proposed  approach,  then,  is  justifiable  in  a 
mathematical  sense.  The  Justification  is  with  respect  to 
simple  criteria  that  are  pertinent  both  to  the  operator  of 
the  system  and  to  the  passengers.  The  mathematical  sense  is 
asymptotic  and  probabilistic,  and  will  be  discussed  shortly. 

From  a practical  viewpoint,  the  design  approach  has 
many  attractive  features.  Most  important  is  that  the 
schemes  have  a decomposed,  modular  nature.  This  yields  a 
remarkable  simplicity  at  many  levels.  Computational 
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requirenents  are  kept  low:  no  long  searches  are  needed,  and 
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dynamic  behaviour  of  each  bus,  and  hence  to  measure  the 
performance  of  a suggested  scheme  In  terms  of  the  parameters 
(e.g.  the  number  of  buses,  the  size  of  the  region,  the 
number  of  transfer  points,  and  so  on).  One  can  thus  choose 
parametric  values  'optimally'  and  can  evaluate  and  compare 
alternative  schemes. 

As  remarked,  we  have  concerned  ourselves  with  the 
fundamentals  of  algorithmic  design.  Therefore,  we  have  had 
to  avoid  operational  details  - there  are  very  many  practical 
aspects  which  we  have  barely,  if  at  all,  touched  upon.  So, 
there  is  scope  for  much  additional  work.  Further  analytic 
work  could  focus  upon  a specific  methodology  for  choosing  a 
decomposition  of  the  region  and  a specialization  of  the 
buses  in  subregions,  particularly  when  the  demand  is 
noDuniform.  We  have  but  indicated  (in  Chapter  ZV)  how  the 
design  tool  might  be  sharpened  and  calibrated  if  stochastic 
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approach  to  coBbina torial  optimization  suggested  by  Karp 
[9].  Karp's  approach  has  permitted  a global  view  of  the 
problem:  when  a new  demand  arises  there  is  no  need  for  a 
large  combinatorial  manipulation  of  the  existing  solution, 
and  a small  adjustment  suffices.  It  yields  simple, 
decentralized  algorithms  that  are  efficient  (in  terms  of 
computational  requirements)  and  yet  are  nontrivial. 
Optimality  can  be  essentially  guaranteed  when  the  problem 
size  is  very  large.  When  the  problem  size  is  small,  the 
theory  la  no  longer  valid,  but  the  resulting  algorithms  are 
nonetheless  relevant.  Thus,  we  are  led  to  claim  that  Karp's 
approach  does  appear  to  have  practical  worth. 


In  applying  [9]  to  our  transportation 
generalized  a result  of  Beardwood  et. 
tempting  to  conjecture  that  similar 


problem  we  have 
al.  [2].  It  is 
generalizations. 


-104- 


yieldlng  interesting  practical  algorithms,  might  be  possible 
for  other  travelling  salesman-llke  problems  that  abound  In 
Operations  Research. 

Traditionally,  scheduling  has  been  regarded  as  a 
combinatorial  problem.  This  emphasis  has  produced  largely 
Intractable  problem  formulations  and,  except  for  some 
special  cases  (see  Conway  et.  al.  [6]),  few  Insights.  Our 
study  has  essentially  avoided  combinatorics.  It  Is  believed 
that  many  of  the  features  we  have  obtained  are  more 
generally  relevant  to  scheduling  applications.  Particularly 
when  there  are  stochastic  disturbances  (as  Is  the  case  In 
most  large  real  problems)  It  Is  felt  that  combinatorial 
aspects  can  be  dissipated  In  a fairly  simple  way,  by 
decomposing  and  discretizing  the  problem  suitably.  Then, 
other  effects  predominate:  these  are  problems  of 
coordination  and  of  designing  the  schedule  In  advance.  The 
schedule  should  be  'stable'  with  respect  to  reasonable 
disturbances  and  should  be  'flexible'  so  that  changes  can 
easily  be  made  when  required.  Phenomena  such  as  these  are 
not  well  understood  within  the  present  framework  of 
Scheduling  Theory. 


AfPEMCIX  CHAPTER  II 


This  appendix  is  a collection  of  detailed  proofs  of  the 
assertions  In  Chapter  II.  All  notation  and  terminology  Is 
as  described  there.  For  the  most  part  the  results  are 
proved  rigorously;  for  the  sake  of  clarity  in  Proposition 
3.2  we  do  revert  to  a more  informal  presentation. 

First  we  make  some  preliminary  notational  remarks,  and 
restate  the  law  of  large  numbers  and  Beardwood's  result  for 
completeness . 


fftliBlniry  aanArka 

Throughout  the  appendix  we  shall  be  dealing  with  random 
variables.  We  must  explicitly  identify  the  underlying 
probability  space  and  describe  our  notation. 

Let  w s { ( o^ , d ^ } , ( 02 1 d2 ) • • ‘ * countable  sequence 
of  pairs  of  points  in  the  planar  region  R.  Let 

''n  • ( (oi  ,di ) , . . . (on»dn)  J be  the  first  n pairs  of  w. 
The  set  W of  all  sequences  w is  the  sample  space.  If  we 
regard  w as  a sequence  of  'random'  pairs,  each  element  of 
{o^ , d ) , 02 t d2 I 03 , d3 , . . . ) being  independently  and  uniformly 
drawn  from  the  points  in  R,  then  we  have  a probability 
space.  Random  variables  are  all  defined  on  this  space.  The 
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set  of  'problem  instances'  of  size  n is  {wn»  w€W}. 


For  example,  is  the  length  of  the  optimal  bus 

tour  through  the  first  n points  of  wCW,  and  is  a 

random  variable.  If  Xjj(‘),  n=1,2,...  are  random  variables 

then  limsup  Xn(')  and  liminf  X^C")  are  well  defined 

n->aD  n->OD 

random  variables.  If  X and  Y are  random  variables, 

X^Y  means  X(w)iY(w)  for  all  wOW;  X^Y  almost 
everywhere  la.e.)  means  Pr L X ( w)iY ( w ) ] = 1. 

For  clarity,  in  what  follows,  we  shall  omil  explicit 
mention  of  the  sample  space  and  the  parameterization  on  the 
w' s. 

Kolmogorov ' s StfOPK  L&H  Sil  Large  tjuflJBACa  f 14 1 


Let  {Xj^,  lsl,2,...}  be  a sequence  of  Independent 

Identically  distributed  random  variables,  such  that 

E (X  1 ) <00.  Then 

n 

J.  51  X^  converges  almost  everywhere 


lim  J. 
n-^oo  n 


n 


5:  Xi  = 
in 


E(Xi  ) 


a . e . 
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2^  Hu  Travelling  Salesman  £roblem 

Theorem  1 (Beardwood  et.  al.  [2]) 

Let  R be  a planar  Lebesgue  set  with  area  a>0  and 

let  Ln  be  the  length  of  the  shortest  path  through  n 

points  which  are  uniformly  and  Independently  distributed  in 

R.  There  exists  an  absolute  constant  b,  Independent  of 

the  problem  instance  and  of  the  shape  of  R such  that 

11m  ^ t bTa  almost  everywhere  (a.e.) 

n^oo  Tn 

Let  us  prove  an  easy  Initial  lemma  that  will  be 
required  in  many  proofs  to  follow:  if  the  number  of  points 
in  a certain  subregion  is  random,  we  can  replace  this  random 
number  by  its  asymptotic  value  in  Beardwood's  formula. 

Lemma  1 . Q 

Consider  a particular  region  of  area  a in  which  there 
Pn  points;  Pn  is  a random  variable  satisfying  for 
some  constant  g, 

11m  £iiEg  a.e.  (Al.O) 

n->^0D  n 

Let  hi)  be  the  length  of  an  optimal  travelling  salesman 

tour  on  these  Pq  points.  Then, 

11m  s byg^a  a.e. 

n->0D  /n 

ErooX 

By  Theorem  1,  we  have 

iim  jLn_  » pyi  a.e.  (A1.1) 

Pd^oo  yp^j 
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By  (A  1.0), 


11m 

n-^OD 


PJI 

n 


Vg 


a . e . 


Let  any  c>0  be  given.  Let  c 


1 


(A1 .2) 


max { e /3 * ^ /BbTa ) and 


Eg  * max{e /3 ,e /37g} . From  (A1.2)  we  have  that  there  exlsta 
Ki  such  that  for  n2N i , 

i i (yg+e^)yn  a.e. 

From  (A1.1),  there  exists  N2  such  that  for  n2N2, 

(bya-e2)y'i^  ^ hn  < ( by^^  ^ a.e. 

Thus,  for  n 2 max{Ni,N2}, 

/ . X s V 

e • 


(bya- E XTg-  ) < iJji  < (byi+e  )(yg+c  ) a. 

^ ’ yF  2 1 


1 . e . 


byayg  - e < iln 
yn 


byiayg 


+ e 


a.e. 


[] 


and  the  result  follows. 


2.2  Static  aiD£l,fi.=fiUS  Problem 


QbaerYatica  2...fl 

The  single-bus  static  problem  is  NP-complete. 

Proof  (See  Aho  [1]  for  the  notation.) 

It  Is  easy  to  see  that  the  bus  problem  Is  in  NP. 

Given  an  instance  of  the  TSP  of  size  n,  we  can 
transform  it  In  polynomial  time  into  an  instance  of  the  bus 
problem  of  size  n as  follows.  Let  the  given  n points  be 
origins  {01,... Oq},  and  let  all  have  as  common  destination 
s point  p,  a distance  of  at  least  from  every  origin, 

where  0 is  a very  large  real  number  and  A is  the 


From  Figure  A2 . 1 , 
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2m  2m  2m  

f ya  dx  <.  zyi  s.  / yx+i  dx 

0 i=l  0 


i . e . 


i . e . 


1 . e . 


m 


^71 


2(2m) 

3 


3/2- 


2m 

I yr  i 


.3/2  i = 1 


l(2m+1 


£23/2  ^ 


2m 


I.  yi  1 £(2  + l/m)3/2  _ 
3/2  1=1  3 

2m 


3a 


3/2 


iy?  i 1 I yr  i ly?  ♦ od/m) 

3 „3/2  i=1  3 


[] 


Lemma  2.2 


limsup  La  i Jiy^pya  a.e. 

n-^os  yii  3 


Proof 


Conaider  the  feaaible  bua  tour  given  by  Algorithm 


Let  Oj^  be  the  number  of  origlna  in  each  region  rj^ 

d^  be  the  number  of  deatinatlona  in  rj^  that  are 
viaited  on  the  firat  paaaage  (with  orlgina  in 
reglona  r i , rg , . . . r^. i ) 

S 

dj^  be  the  number  of  deatinatlona  in  rj^  that  are 

viaited  on  the  aecond  paaaage  through  the  reglona. 
Let  a^  -be  the  length  of  the  ahorteat  tour  through  the 
o^-fSi  polnta  in  rj.  viaited  on  the  firat  paaaage,  and 

a^  be  the  length  of  the  ahorteat  tour  through  the  d^ 
pointa  viaited  on  the  aecond  paaaage  - theae  are  random 


varlablea . 
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If  Yjj  is  the  length  of  the  optimal  bus  tour  then 
m . m 


“ 1 “ ? 

Yn  Z Si  Z Si  ♦ 2af^(R} 

i=1  i=1 


(A2. ^ } 
(so  A(R  ) /Jn  -»  0 as 


where  A(R)  is  the  diameter  of  R 
n-»0D) . 

Now,  the  law  of  large  numbers  yields,  as  n-^oo^ 

a . e . 


£1 

n 


and 


1 

o 

ii-A) 

2 

m 

ffl-l.tl 

2 


Thus,  by  Theorem  1 and  Lemma  1.0, 

b 


1 2 

lim 

n->®  Vn 


1 -■ 

m 


a . e 


a . e . 


A 4-  b 

m 


a . 6 


* ( ^m*  1 ♦ Vm*  1^1) 

_3/2 


a . e . 


Hence, 

" 1 

lim  _L  [ Z aj 
n-^®  ^n  i»  1 


“ p 

♦ I .f) 

i«l 


2m 

= -bi/a  [ z yr  y*  - yjm] 

-3/2  1»1 


a . e . 


Inequality  (A2.1)  now  yields 

2m 

limaup  la  i .b^a  ( Z yi)  ♦ b( i-y?!  yi 

n-^®  yii  _3/2  i«1  a 

■ 

This  holds  for  arbitrary  m,  so  we  can  let  m->oo»  and  we 
obtain  by  Lemma  2.1, 

a . e. 


limaup  la  1 lyjbyi 
n-^®  yij  3 


[] 
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llfflinf  JLq  2.  Jl>/?bya  a.e. 

n-^oj  yn  3 

IXQSil 

First,  some  notation.  In  any  problem  Instance,  let 
On  be  the  set  of  n origins  and  Dn  be  the  set  of  n 
destinations.  Let  Tn  be  any  optimal  tour  (which  has 
length  Yn).  By  this  we  mean  that  Tn  is  that  connected 
subset  of  R which  is  traversed  in  visiting  On  U Dn*  A 
tour  segment  through  P (where  P c Tn)  is  the  minimal 
connected  subset  of  Tn  that  contains  P. 

The  optlm>al  tour  Tn  induces  a feasible  ordering  on 
the  points  of  On  U Dn-  Let  us  label  the  origins  in  On 
such  that 

KJ  <s>  ©i  proceeds  oj  on  the  tour  Tn- 
Let  d^  be  the  destination  corresponding  to  o . 

Let  m be  any  Integer,  m>1.  We  partition  R into 
m Lebeague  subregions,  each  of  area  a/m. 

Let  (Yn)^  be  the  tour  segment  through 
{oi ,02, . . . } , where  [n/m]  is  the  smallest  integer 
larger  than  n/m.  (The  segment  (Tn)^  possibly  visits  some 
destinations  en-route.)  Similarly,  let  (Tn)'^  be  the  tour 

segment  through  ^ «>[  ( j-i  )n/m]*1 ' * * * ®[  Jn/m]  ^ 
partitions  of  R into  m Lebeague  subregions  Ri,...R| 

eaoh  with  area  a/m  and  bounded  perimeter  (length  of 
boundary)  such  that  the  origins  in  Rj  are  precisely  those 
in  (Tq)^,  for  jBl,...m.  Let  A^,...Ab  form  a partition 
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of  R drawn  at  random  from  the  space  of  such  partitions ^ ^ ^ . 
We  can  assume  that  the  origins  in  Aj  are  uniformly 
distributed  there. 

Let  us  define  '^r  each  subregion  Aj,  Jsi,2,...m, 

Onj  s number  . rigina  in  Aj  x [n/m]  or  tn/ml+l 

dnj  s number  of  destinations  in  Aj 

d^j  * number  of  destinations  in  Aj  which  also  belong 

to  ^ 1 • **2 »•••<*[  jjj /m]  ^ " i.e.  destinations  in 

Aj  whose  corresponding  origins  are  in 
Ai  U A2  U ...  U Aj. 

(1)  An  example  of  suitable  Aj , Jx1,...m,  is  provided  by 
the  following. 

Let  ix1,...n,  be  a ball  around  Oj;^  with  radius 

c<l/n^  and  with  £1  Bj  x 0 for  i^J.  That  is, 
c < minll/n^,  mini lloi-ojll/2 , i,Jri,...n,  ii^J)} 

and  Bi  X R n lx,  llx>Oill<e),  ixl,...n. 

m 

Now  partition  R-^U^Bi  randomly  into  a equal  regions 
such  that  each  Xj  is  a Lebesgue  set  with 
perimeter  leas  than  a fixed,  suitably  large  number,  K. 

Finally,  let  Aj  be  Xj  together  with  the  balls 
around  the  origins  in  (Td)^,  i.e. 

*J  * ®[(J-1)n/m34.i  •••  “ ®[Jn/m] 

The  perimeter  of  Aj  is  less  than 

2 

I ♦ n.2v/n  * * ♦ 2v/n 

(This  is  a bounded  function  of  n,  as  is  to  be  required  in 
(A2.6a).) 
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**  n 1 

d^j  - number  of  destinations  in  Aj  fl  (T^)  * l.e. 

destinations  which  are  visited  by  the  segment 

(Tn)^  in  Aj. 

Note  that  d^j  ^ <inj' 

Since  the  regions  Aj  were  randomly  selected,  and  the 
destinations  were  uniformly  distributed  over  R,  the  random 
variables  d^j  and  d^j  satisfy,  by  the  law  of  large 


numbers , 


Also , 


n 

Ik 

n 

®nj 


J.  a . e . as  n^oo 

m 

^ a.e.  as  n->0D 

m2 

1 as  n->OD 

m 


such  that, 

(A2.2) 

(A2.3) 

(A2.iO 


Now,  the  tour  within  Aj,  viz.  Tq  n Aj,  consists  of 
one  or  more  connected  pieces.  Discard  any  such  piece  which 
does  not  contain  a point  of  Oq  U D^.  The  closure  of  any 
surviving  piece  will  be  called  a J-piece.  The  set  of 
j-pleoes  can  be  partitioned  into  two  parts  (see  Figure 
12.2): 

^•pieces  are  those  that  belong  to  (T^)*^ 

(2) 

J -pieces  are  the  remainder. 
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Figure  A2.2  The  Tour  T*-(i) 

Port  of  the  optimal  tour  T*  is  shown.  The  arrows  indicate  the 
direction  traversed. 

j^’^-pieces  ore  shown 

/2)_  pieces  are  shown  + + + + + + 

(Tp)^  ” segment  is  the ond  ==  path. 


Thus 

, the 

-pieces  contain  all 

Onj  origins  In 

Aj  as 

well 

as 

the 

d^j  destinations. 

The  J^^^-pleces 

contain 

dnj 

- dnj 

destinations . 

Let 

be  the  length  of  the 

J ^ ' ^-pieces 

.12) 

be  the  length  of  the 

1 2 ) 

j'  -pieces. 

Clearly, 

In  . 

i 1 (sj^^  ♦ sj^^ 

J-1 

(A2.5) 

Let 

.•(1) 

be  the  length  of 

an  optisal  travelling 

•slesaan 

tour 

through  the  Onj 

* <^nj  points 

of  the 
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,( 1 ) , 

J -pieces. 

M ( 2 ) 

Let  sj  be  the  length  of  an  optimal  travelling 

salesman  tour  through  the  d^j  - ^nj  points  of  the 

(2) 

j'  -pieces. 

For  k=1,2,  Beardwood  et.  al.  (their  Lemma  2)  have 


shown  that  there  exists  a Aj  such  that 


as  n-^oo 


(A2.6a} 


satisfying 


• (k)  . _(k)  y.(k) 

sj  1 sj  ♦ Aj 


(A2.6) 


This  is  obtained  by  constructing  a feasible  travelling 

(k) 

salesman  tour  through  the  points  of  the  J -pieces,  each 
traversed  Just  once,  together  with  a part  of  the  double 
circuit  of  the  boundary  of  Aj.  (An  example  is  diagrammed 
in  Figure  A2 . 3 • ) 

Now,  there  exists  an  N2  such  that,  for  all  n2N2« 

Ak) 

0 < Zj < £ for  k*1,2;  J=1,2,...m.  (A2.7) 

■/E  m 

By  Theorem  1,  for  n>N2, 


b A - i i 
Jm  m 


.•(  1 ) 

"J 

^Onj  ♦ 


(A2.8) 


bfi  - i i 


m m 


.•(2) 



<5nj  - 


(A2.9) 


From  (A2.8),  (A2.9),  (A2.2)  and  (A2.<l),  we  get  for 

d>N  « mmx(N^ ,N2) , 


Figure  A2.3  The  Tour  Tn  -(ii) 

(2) 

j -pieces  ore  shown  +++  + -♦-  + 

The  tour is  on  upper  bound  to  the  length  s**^\  on 

optimal  trovelling  salesman  tour  through  the  destinations 
in  the  pieces. 


It  is  true  in  general  that,  for  any  a>b>y>0, 


argain  '4>  yb-x}  « y 

xe[0,y] 

So,  ainoe  dQj  e [O.d^j],  and  vltb  (A2.3)  we  have 


(k2.9a) 


But,  e is  arbitrary,  and  the  coefficients  of  e ar« 
bounded  above,  so 

n 

liBinf  In  1 b^a  I (yE+3  ♦ v^5^)  a.e. 
n->0D  Vn  -3/2  Jsl 

D 


This  holds  for  any  m.  Thus,  by  Lemma  2.1, 
liminf  2Ln  1. 


a.e. 
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^.tatig  MultlPle-Bua  Problema 


Lgmna  1 

-k 

lim  ±ii  - cTa  a . e . 

n->oo  yH 


Proof  _ 

Define  the  random  variablea 

-k  7^ 

e = limaup  tsi 
n->®  Jn 

k 7^ 

£ = llmlnf  ia 

n->cD  Vn 


and 


k 


£. 


aa  followa. 


Note  flrat  that  It  la  feaalble  to  have  k-1  buaea  Idle 
and  to  uae  only  one.  Thua , 

Zn  ^ Yn 

(Yf]  la  aa  In  Theorem  2.) 

a.  k 

Hence,  c ^ c a.e.  by  Theorem  2. 

If 

We  now  ahow  that  c a.  a.e. 

With  the  fact  that  i c**,  we  ahall  have  the  required 
reault. 

Define  for  ls0,1,2,...  auxllliary  problema  aa 

followa . 

^Problem  Let  ^j^d)  be  the  minimal  total  dlatance 

travelled  by  k buaea  through  n random  demand  paira, 
allowing  "Qat  1 tranafera  of  paaaengera  - Y|j(l)  la  a 
random  variable. 

Let  limaup  ^n^ ^ ^ 

ar^CD  n/a 

liminf  ) 
n^oo  yn 

Now,  for  la0,1,2,...  , we  have 


* 5‘'(i) 
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Y^d)  2 Yn(l  + 1) 

and 

Y^li)  2 Yg. 


So , 

and  also 

£*'(1)  2 S.^  for  all  i. 

Thus,  the  (i  (i),  i=0,1,2,...}  form  a monotonically 

decreasing  sequence  in  i,  bounded^J)elow  by  £ . Thus, 

lim  i'^Ci)  = exists,  with  2 

i->(r) 

It  is  shown  in  the  Lemma  3.1.1  below  that 

(A2.10) 


For  any  fixed  i,  consider  the  optimal  tour-length 
Yji(i)  - a random  variable.  It  is  possible  to  Justify  - 
see  Lemma  3.1.2  below  - that 

Yn(i)  ♦ (i+1)(k-1)A  2 yJ  (A2.11) 

where  Y^  is  the  optimal  distance  travelled  by  a single 
bus,  and  A is  the  diameter  of  the  region  R. 

Dividing  (A2.11)  by  and  letting  n->ac^  we  obtain 


by  Theorem  2 that 
ll‘‘(i) 

This  holds  for  any 

2 

and,  by  (A2.10), 


2 c a . e . 

i;  thus, 


c a . e . 


c a . e . 


This  proves  the  result,  modulo  the  two  lemmas. 
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L£mm&  3.1.1 

T k k 

Vie  need  to  show  that  { = £.  • Suppose  to  the  contrary 


that  > XL*'  at 

some  wBW.  Let  c : 

- > 0.  Now, 

there  exists  an  1 

such  that  for  i2I» 

J*'  - e/6 

< U ) 

Also,  since 

liminf  7nli)^ 

n->oo  ^n 

there  exists  N such  that  for  all  n2N, 

1 - £ 

6 

- see  Royden  £20,  page  37]. 

Thus,  for  121  and 

n2N, 

1 - 1 

yh 

3 

= ♦ 2 e 

3 

(A2. 12) 

Since 

7k 

liminf  hi 

k 

' s.  . 

n->0D  /n 

there  exists  an  h>N 

such  that 

in  < il 

♦ e 

(A2. 13) 

yR 

3 

- again,  see  Royden 

[20,  page  37]. 

But  Zb  is  the 

length  of  an  optimal 

tour  on  h 

demands.  This  tour  can  have  only  a finite  number  of 

transfers,  so  there 

exists  an  integer  Ph^c) 

suoh  that  for 

almost  every  problem 

Instance,  z{j  requires 

at  most  Pb(^) 

transfers.  Thus,  for  i>Pb(^)» 

7k 

Zh  * 

»h(i) 

(A2. 14) 
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Kence,  for  n2N  and  i 2 maxll,  pjjCc)}.  we  have  from 

(A2. 12)  - (A2. 14)  , 

V If 

i ♦2^  < £ ♦£. 

3 3 

This  Is  clearly  absurd  and  the  lemma  is  proved. 

Lemma  5>.K2 

We  assume  that  the  buses  start  their  tours  from  certain 
fixed  (but  arbitrary)  points  in  R,  and  that  all  buses 
terminate  their  tours  at  a common  final  point  in  R. 

Consider  the  following  acyclic  directed  graph,  G, 

representing  the  tour  with  length  problem 

instance . 

Represent  by  nodes  the  k starting  points,  the 

transfer  points  (there  are  at  most  1 of  them)  and  the 
final  terminal  point.  The  tour  of  a particular  bus  i is 
represented  by  a path  through  G:  edges  into  a node 

represent  the  buses  which  transfer  passengers  there,  and  the 
indegree  equals  the  outdegree  at  all  nodes  representing 
transfers.  An  edge  thus  represents  the  path  through  R 

taken  by  the  corresponding  bus,  and  G represents  the  given 

If 

set  of  tours  with  total  length  Y^d). 

An  example  best  illustrates  this  definition  of  G.  The 
graph  In  Figure  A2.4  represents  a tou.*'  with  7 buses  and  6 
transfer  points. 

Suppose  the  k tours  are  executed  by  a single  bus  In 
the  following  way.  Pass  along  the  same  routes, 
' backtracking  In  a straight  line  from  each  node  to  a 
previous  one  as  far  as  Is  necessary,  so  that  all  Incoming 
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edges  have  been  executed'.  (We  do  not  feel  that  it  Is 
necessary  to  be  more  specific  here.  The  idea  is  precisely 
the  reverse  of  the  standard  depth-first-search  algorithm 
see  e.g.  Aho  et.  al.  [1].)  At  worst  we  will  have  to 
backtrack  k-1  times  at  each  of  i-«-1  nodes  (the  transfer 
points  and  the  terminal  point),  and  each  backtrack  will 
involve  a distance  of  at  most  A,  the  diameter  of  R. 

This  leads  to  a feasible  single-bus  tour  with  length 
» S l5(l)  ♦ (k-lXl  + DA 


[] 


yielding  inequality  (A2.11). 
Lemma  3.1  Is  thus  proven 
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Theorem  2 

If  transfers  are  allowed 

lim  in  = a.e. 

n->oo  y/n  k 

Proof 

i 

£ Xn  Is  the  total  distance  travelled  by  the  optimal  tour. 
1 = 1 

ir 

If  Zn  Is  as  In  Lemma  3-1f  then 

If  1 k 

^ *n  ^n 

1 = 1 

. k 

Now,  kin  i Z Xn  by  definition  of  Y^. 

i=l 


Dividing  by  Jn  and  letting  n-^oo,  we  obtain,  by  Lemma 

3.1, 

1 

limlnf  *^^n  2 cTa  a.e.  (A2.15) 

n-^oo  yn 

k® 

Consider  the  tour  given  by  Algorithm  3.  For  this 

tour  let 

Sn(m)  be  the  total  distance  travelled  by  bus  i 
hn(m}  be  the  total  length  of  the  travelling  salesman 
tours  covered  by  bus  i in  its  regions. 

Then, 

*i(m)  i 

where  A is  the  diameter  of  the  region  R. 

Now,  by  Lemmas  1.0  and  2.2, 


llmsup 

hi(£) 

S.  mcflfn  * 0(1/m)  a.e. 

n.>a) 

yn 

|mja 

Hence, 

llmsup 

^ cyi  ♦ Od/m)  a.e. 

n^oo 

yn 
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By  the  optimality  of  Yn  have 

Yn  i max{sn(m),  is1,2,...k}  for  any  m. 

Thus , 

limsup  1 c7a  ♦ 0(1/m)  a.e.  for  any  m. 

n->OD  Jn 

Letting  m->ao, 

limsup  in  i a.e.  (A2.16) 

n->OD  ^n 

The  inequalities  (A2.15)  and  (A2.16)  prove  the  lemma. 

[] 

ProDOaition 

The  3>bus  fixed-route  scheme,  described  in  Section 

3F 

2.3.3  yields  a time-to-delivery  Yj;  satisfying 
3^ 

11m  ^n  s 1 (yj  ♦ yj  ♦ 1)bya  a.e. 
n->oo  yiT  3yT 

• .798bya  a.e. 

Proof 

Ue  use  the  same  notation  as  that  used  for  describing 

the  algorithm  in  the  text.  Our  derivation  will  be  Informal, 

but  hopefully  transparent.  We  need  onlv  consider  bus  1 

since  the  other  buses  have  analogous  tours  and  so  distances 

travelled  will  be  the  same  (asymptotically  a.e.). 

Bus  I first  visits  regions  ( r ^ , r2  • . . . rji ) . Zn  each  r^ 

2 

[ (a-i'«>1  )4>a]n/9a  origins  are  collected  and 
2 

[i-l]n/9a  destinations  are  delivered. 

80,  by  Lemma  1.0,  there  is  a tour-length  of  bya/3my^n/9m 
in  ri . 
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Second,  the  bus  visits  ( r2D • r2B- 1 • • • • ‘'■•f  1 ) • each 

2 

[ m*  ( 1 - 1 ) ] h/Sjb.  origins  are  collected  and 

2 

[ (m-1  )-*-m] n/9ni  destinations  are  delivered. 

The  tour-length  In  each  of  these  m regions  is 
bVa/SniJ  ( 3ni  - 1 ) / 9ni^yn . 

Finally,  the  b\ia  visits  (ri,...r]|))  again.  In  each  r^^ 
2 

(a-«-l}n/9ni  destinations  are  delivered. 

So  in  r^  there  is  a tour-length  of  bya/SnJ  (m  + 1 )/9ni^yn. 

The  total  distance  travelled  by  the  bus  is  thus 

byjyniyi  yr-TTm  ♦ yr+TTE) 

3^1 

* byaTnCy?  + ♦ 1 ) •*•  0(1/in)yn 

3^1 


[] 


Letting 


m-^oD  we  obtain  the  required  result. 


APPENDIX  CHAflEB  III 


This  appendix  collects  proofs  of  the  assertions  in 
Chapter  III.  In  contrast  with  Appendix  II,  proofs  are 
inforaal.  With  additional  work  the  results  can  be  made 
rigorous  • for  example,  Lemma  2.2  can  be  used  as  a model  for 
Lemma  4.1  and  Lemma  2.3  as  a model  for  Theorem  4.  However, 
it  is  believed  that  this  would  blur  their  simplicity. 


3. 1 The  Static.  Sinale-Bus  ProblCB 


Lemma  4 . 1 


(1)  lim  Wn(Ta) 
n->OD  >/n 


(ii)  11a  ^n^T,) 
n->co  Vri 


2 [8y?  - 7]  bVa  ♦ 0(1/m)  a.e. 

15 

.5752bya  ♦ 0(1/a)  a.e. 

1 [118^7  - 47]  bTa  ♦ 0(1/m)  a.e, 

105 

1 . 1417bya  Od/m)  a.e. 


graaX 


Label  the  a regions  ri,r2»...ra. 

By  Leaaa  1.0,  we  can  assuae  that,  for  any  l,js1,...m. 


are 

n/a^ 

deaanda  of  the 

fora 

(ri,rj).  Let 

ti 

be  the 

tiae  spent  in 

Ti  on 

the  first  passage 

ti 

be  the 

tiae  spent  in 

Ti  on 

the  second  passage 

I ■ 

a 

^ ti 

> tlae>to>ooaplctlon 

of  first  passage. 

Ia1 
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For  n large  enough  we  have,  with  probability  1, 

t i * b^n  ym+1-1  ( A 3 . 1 ) 

D 

2 

since  n/m  origins  and  (l-1)n/m  destinations  are 
collected.  We  can  assume  that  the  expected  time  at  which 
any  of  these  points  is  visited  is  ( 1 1 +t2+ • • . ••■ti_ i ♦t  /2  ) . 
Similarly , 

ti  = b^n  ym-i+1  (A3. 2) 

m3/2 

2 

since  (m<-i  + 1)n/m  destinations  are  visited  in  each  region 

r^.  The  expected  time  at  which  these  points  are  visited  is 
T 4-  ti  + t2  ...  t[/2. 

First,  let  us  consider  the  waiting  times.  The  sum  of 
the  waiting  times  in  region  1 is 
i-1 

H [ Z tj  ♦ ti/2]. 
m J.l 

The  total  average  waiting  time  is  then 

m 1-1 

w^CTb)  X JL  ja  Z [ Z tj  ♦ ti/2] 
n m 1x1  Jxi 

m 

X byii  [ Z (m-i4-1  / 2 )ym4-i-1  ] 

7/2  1x1 

D 

by  expanding  and  using  (A3.1}. 

Using  a technique  similar  to  that  of  Lemma  2.1,  it  can 
be  shown  that  this  series  converges  and 

■ JL  [8y?  - 7]  ♦ 0l1/m) 

15 

= .5752  ♦ 0(1/m) 

Equation  (ii)  is  somewhat 


1 Z (a-i4>i  /2 )ya-*'i.i 
.3/2  1x1 


This  proves  Equation  (i). 
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I 

I 

i 


■or«  cumbersome. 

For  destinations  which  are  visited  on  the  first 
passage,  the  sum  of  delivery  times  Is 
m 

G s Z (i-l)ja_  (ti  ♦ ...  ♦ ti.i  ♦ ti/2) 

’ 1*1  2 

B 

* J,  Z ym-4-i-l  (m^  - m-  1^  + 21-  1) 

7/2  2 1*1 

B 

Again  it  can  be  shown  that 

°l  * b ri(2^''^-1)  - 1(2^'^^-1)“|  ♦ 0(1/m) 

„3/2  [5  7 J 

The  sum  of  delivery  times  for  passengers  delivered  on 

the  second  passage  is 


m 

G«  a Z [T  t|  ♦ ...  ♦ tj_^  ♦ tl/2}  (m-l-»1  )n 

^ i«1  .2 

* bn^*^^  [ Z ym+i-1  ^(■'*'1)  ♦ ^ (m-i<f1)^^^] 

^7/2  1*1  2 1*1 

Once  again  it  can  be  shown  that 

°2  * b ri(2^^^-1)  ♦ ll  ♦ 0(1/m) 

„3/2  [3  7j 

Finally,  the  average  flow  time  is  (G^<«-G2)/n,  l.e. 

* b r^(2^^^-1)  - 1(2^^^-1)  ♦ 1(2^^^-1)  ♦ ll  ♦ 0(1/m) 

[5  7 3 7j 

* _h_  [IlSy?  - A7]  ♦ 0(1/m) 

105 

= 1.1417b  ♦ 0(1/m).  [] 

The  proof  of  Lemma  4.2  follows  easily  from  Lemma  4.1  { 

and  the  dlsousaion  in  the  text  of  Section  3>1«  { 


3 


IheoceB  Ji 

Let  fj  = inf{fn(T),  TCS}. 

Given  any  OO,  there  exists  an  M(c)  such  that  for 
iii2M  and  for  n large  enough,  la  e-optlmal  for 

prohlea  ( 3 • 1 ^ a . e . 

Hence , 

11m  = 1 . UlTb^a  a.e. 

n->oo  Tn 

Lrool 

Throughout  this  proof  we  assume  that  n Is  very  large, 
so  results  obtained  are  asymptotic  In  n.  For  simplicity, 
we  drop  the  parameterization  on  n. 

For  any  problem  Instance,  let  the  tour  Tq  minimize 
{f(T),  TeSj,).  Thus,  for  all  m, 

f(TS)  1 f(TS). 

Let  OO  be  given,  and  suppose  that  the  hypothesis  Is 
not  true.  Then,  for  a subset  w'  of  W with  non-null 


measure,  there  exists  such  that  for  all  m2.M^, 

f(TS)  - f(TS)  > 3e/2  (A3. 3) 

By  Lemma  4.3  there  exists  such  that  for  all  m^N^, 

w(tS)  ♦ fdS)  - w*  * f*  i e/2  a.e.  (AS.'*) 
So,  for  a > aaxiM^,M2}, 

fdj)  ♦ w(tJ)  2.  f*  * V*  a.e. 

2 w(t2)  ♦ fdS)  - E/2  a.e.  by  (A3. 4) 

2 w(tS)  ♦ fdS)  ♦ c a.e.  by  (A3. 3) 

i.e.  wdj)  - wdS)  1 E A.t.  in  W’  (A3. 5) 

For  any  partition  of  R into  m equal  subregions  we 


oan  define  a (nonfeasible)  tour 


that  is  similar  to 
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O 0 

Tji:-  visit  the  regions  successively,  as  for  tour  Tm,  but 
in  addition  visit  the  destinations  of  the  form  (rj^,r^)  in 
each  region  r^  on  the  first  passage  through  the  m 
regions.  Clearly, 

w(sS)  > w(Tb). 

But,  by  choosing  m large  enough  - i.e.  in>M^,  say  - we 

w(sj)  - v(tS)  < (A3. 6) 

Let  M s maxlM.j  . For  any  m^M  consider  the 

tour  Tg  for  any  problem  Instance.  We  construct  another 
partition  of  R.  (This  is  again  an  extension  of  the  method 
in  Theorem  2;  see  that  proof  for  notation.) 

Divide  the  tour  into  two  segments:  S*  is  the 

segment  from  the  start  of  the  tour  to  the  delivery  of  the 
last  origin;  S"  is  the  remaining  segment  containing  only 
destinations.  Now  divide  s'  into  m subsegments: 
is  the  first  segment  from  the  first  origin  to  the  (n/m)th 
origin;  s|  is  the  (l)tb  segment  from  the  ((i-l)n/m)th 
origin  to  the  (in/m)tb  origin.  There  exists  a partition 
of  R into  m subregions  lr^,...rB}  such  that 

n s|  W 0 (r^  is  the  interior  of  r^) 

r®  n sj  » 0 for  iWJ 

(Note  that  s'  might  cross  Itself,  in  which  case  different 
segments  have  points  in  common.) 


(1)  Observe  that  this  is  precisely  the  technique  used  in 
Theorem  2,  Equation  (A2.9a).  The  relationship  (A3. 6)  can  be 
Justified  in  the  same  way;  the  formulas  are  particularly 
cumbersome,  and  would  obscure  the  simplicity  of  the 
argument. 
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The  space  of  all  such  partitions  Is  nonempty.  Let  us 
draw  a partition  uniformly  from  this  space,  and  consider  the 
tours  t2  and  Sq  respectively  on  this  partition  of  R. 
We  can  assume  that  the  origins  and  destinations  are 
uniformly  distributed  within  each  subregion. 

By  construction, 

wdj)  i w(sS). 

So,  by  (A3. 5)  and  (AS.b), 

w(tS)  ♦ e/2  > w(Sn) 

I wdj) 

> w(tS)  ♦ e 

which  holds  a.e.  In  w'.  This  Implies  that  1/2  >1,  a 

contradiction . 

Thus,  the  first  part  of  the  hypothesis  Is  true.  By 
letting  m-^OD  we  obtain,  from  Lemma  4.2, 

11m  * 1.l4l7bya  a.e. 

n-»o>  yE  [ ] 

3.2  The  Dynamic  Slnala-Bua  Problem 

IMMMA  LlA 

If  the  system  Is  in  steady-state,  then  with  probability 
1,  the  tour  T^  given  by  Algorithm  5 minimizes,  among  tours 
in  Sb,  the  average  flow-time  of  all  passengers. 

Proof 

Label  the  demand  pairs  as  (o^.d^),  is  1 , 2 , . . . n, . . . Let 
t^(i)  be  the  tine  of  arrival  of  (oi,di);  let  tCoi)  be 
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the  time  at  which  Oj^  Is  visited;  and  let  t(dj^)  be  the 
time  at  which  d^  is  visited.  The  average  flow-time  of  the 

demands  is  then  

n 

lim  1 I [t(di)  - t li)]  (A3. 7) 

n->a)  n i = 1 

Now  regard  the  system  as  a queuing  system  (see  e.g. 
[11]).  It  is  convenient  to  consider  two  arrival  streams  - 
the  origins  {o^}  and  the  destinations  {dj^}.  The  arrival 
epoch  of  oj  is  t^(i)  and  its  departure  epoch  is  t(oi). 
The  arrival  epoch  of  d^  is  tCo^)  and  its  departure  epoch 
is  t(di).  Thus,  at  any  time  t the  'units'  to  be  'served' 
which  are  present  in  the  system  are  the  points  which  can 
feasibly  be  visited  at  t. 

The  average  flow-time  of  all  units  is 
n n 

P « lim  _L  { ^ [t(oi)  - t.(i)]  1 [t(di)  - t(oi)]) 

n-»oo  2n  i = 1 i = 1 

n 

« lim  Z [t(di)  - t CD]  (A3. 8) 

n-^OD  2n  i*  1 

which  is  half  the  criterion  we  wish  to  minimize  in  (A3. 7). 
So,  minimizing  (A3>8}  will  also  minimize  (A3-7). 

We  seek  a stationary  decision  rule,  X:-  whenever  a bus 
exits  from  a region,  this  rule  will  determine  which  region 
is  to  be  visited  next  and  which  points  are  to  be  visited 
there.  Any  stationary  such  rule  will  Induce  a stationary 
arrival  pattern  on  the  destinations. 

Given  sny  rule  Z,  let 

Nj  ■ sverage  number  of  units  in  the  system 


s average  time  between  arrivals 
Fj^  = average  flow-time  of  units. 

If  the  system  is  in  steady-state,  a well-known  result  of 
Little  [13]  gives  that 

Fx  = a.e.  (A3. 9) 

Note  that  since  X is  stationary,  the  average  time 
between  arrivals  of  units  is  constant,  l.e. 

Tx  = 1/2q  (A3. 10) 

(q  is  the  constant  arrival-rate  of  demands). 

Further,  if  n units  are  visited  in  a region,  then  the 
average  time  spent  on  each  unit  is  h^/n,  where  h^  is  the 
length  of  the  path  visiting  these  n points.  For  n large 
enough,  ^n^'^  will  be  a decreasing  function  of  n with 
probability  1.  Now,  to  minimize  Nx,  the  average  number  of 
units  in  the  system,  visit  those  which  will  take  the 
shortest  time  - i.e.  minimize  h^/n.  Thus,  choose  the 
region  with  most  points  and,  using  an  optimal  travelling 
salesman  tour,  visit  all  points  there.  By  (A3. 9)  and 
(A3. 10)  this  decision  rule  will  also  minimize  Fx  (a.e.). 

With  (A3. 8)  the  result  is  proven.  [] 

IfatorttB  5. 

Let  ^q,m  8e  the  optimal  average  flow-time  for  tours 

in  Sg  when  the  arrival  rate  of  passengers  is  q.  Then, 

liB  ^9 1»  * 2b^a(1+1/B)  a.e. 

q-^OD  q 

Proof 

By  Laaaa  5.1  the  tour  Blniaizes  the  average 
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flow-tiae  in  Sg,,  so  we  calculate  the  average  flow-time  of 
this  tour. 

Let  6 be  the  time  spent  in  each  of  the  m equal 
(2 ) 

areas  . Then,  by  Equation  (3>3), 

lim  fi  s 2b^a.  a.e.  (A3. 11) 

q-^oo  q m 

Consider  a time  interval  [u,6]  during  which  a 
particular  region  is  visited.  During  this  time  qe  origins 
and  q6  destinations  are  visited;  we  can  take  the  average 
time-of-vlsltation  as  t a ©/2.  During  each  of  the  previous 
m Intervals,  of  6 units  each,  qe/m  origins  and  qd/m 
destinations  arrived  at  this  subregion  - here  by  'arrived' 
we  refer  to  the  arrival  of  units  or  points  at  the  queuing 
system  as  described  in  Lemma  5.1.  We  can  take  the  average 
arrival  times  of  these  2q6/m  points  as  t^  s -(l-1/2}6, 
for  ial,2,...m.  Hence,  the  average  flow-time  of  all  points 
is 

1 [e  ♦ 2©  ♦ ...  ♦ m©]  a (m»1 )© 

2qe  a 2 

By  (A3. 8)  of  Leaaa  5.1,  the  average  flow-time  of  passengers 

is  twice  this,  i.e. 

Fq,B  « (a+l)0 

and 

lim  • 2b^a(1+1/B)  a.e.  by  (13.11). 

q->0D  q [ ] 


(2)  In  the  text  of  Section  3.2  we  called  this  quantity  e 

Q • 
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3.»3  Dynamic  Multlnle-Bua  Problema 

Theorem  ^ 

The  tour  described  by  Algorithm  6 is  optimal 

(asymptotically  a.e.)  for  the  k-bus  minimum  flow-time 
problem  when  tours  are  in  and  at  most  h transfers 

are  allowed  per  period. 

|r 

^q;h,m  average  flow-time  resulting  from 

this  tour.  Then, 

pk  2 - 

lim  q •»  h » “ s ba  p + _L(l  + 1/k)  ♦ l(1-1/k)|  a.e.  (A3. 12) 

q-^OD  q [_  mh  ^ J 

Proof 

If 

We  shall  not  prove  the  optimality  of  the  tour  T^.g,.  A 
proof  could  proceed  along  the  lines  of  Lemma  5.1,  and  the 
same  principle  would  still  be  valid:  each  time  that  a bus 
exits  from  a region,  visit  next  the  region  with  most 
feasible  points,  and  visit  all  feasible  points  there.  For 
q large  this  will  yield  the  tour  of  Algorithm  6.2. 
Further,  with  probability  1,  each  bus  will  spend  the  same 
time  in  each  region,  and  there  will  be  no  waiting  at  the 
transfer  point. 

Below,  we  derive  (A3>12)  in  detail. 

Each  bus  visits  m regions  between  any  two  transfers. 
Let  6 be  the  time  spent  by  a bus  In  any  subregion.  The 
period  Is  then  hm^.  During  this  time  there  are  qhmO  new 
dsaands,  and  so  each  bus  Bust  serve  2qbae/k  points.  In 
aaob  region  there  are  EqhaO/kha  ■ 2q0/k  points  visited. 
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Hence , 


1 . e . 


e 


m? 

k 


2b^qa 

k^ha 


kho 


(A3.  13) 


Each  time  a region  Is  visited  q6/k  origins  and  q6/k 
destinations  are  visited. 

Consider  a bus  B visiting  a region  that  is  the  J-th 

2 

after  a transfer  point.  Of  the  qO/k  destinations,  qO/k 

were  collected  (l.e.  their  corresponding  origins  were 

2 

visited)  by  bus  B itself,  and  (k-1)q6/k  were  collected 

2 

by  the  (k-1)  other  buses.  The  q6/k  points  were  visited 
uniformly  during  the  last  hm  time  intervals,  of  length  0 
each.  So  the  waiting-time  for  these  points  is 

(0  20  ♦ ...  ♦ hmO) 


k^bm 

« _a®.  (hm*1  )0 
2k^ 


(A3. 1«) 


The  (k-1)q6/k  points  were  collected  uniformly  during  the 
hm  time  intervals  before  the  last  transfer,  (J-1) 
intervals  previously.  Thus,  the  total  waiting-time  for 

these  points  is 

>0^  [JO  ♦ (j+1)e  ♦ (J+2)0  ♦ ...  ♦ (hm^J-1)6] 
k^hm 


>0^  (hm+2J-1)0 

.2 


(A3. 15) 


The  qO/k  origins  all  arrived  during  the  last  hm  time 

intervals.  The  total  waiting-time  for  these  points  is 

_3L  (0  ♦ 20  ♦ ...  4.  hmO) 
khm 


• (hmt-l  )e 
2k 


(A3. 16) 
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The  total  waiting-time  for  all  points  visited  during 

this  J-th  region  of  bus  B is  the  sum  of  (A3-14).  (A3. 15) 

and  (A3. 16).  This  is 

_EL.©  [(hm+l)(k+l)  ♦ (hm+2J-1 ) (k-1  ) ] 

2k^ 

and  2q6/k  points  were  visited. 

Summing  now  for  Jsl,2,...a  and  dividing  by  2qem/k 
we  obtain  the  average  waiting-time  of  all  points,  viz. 


a 

_fi_  [m(ha+1  ) (k+1  ) ♦ I ( hm+2 J-1  ) ( k-1  ) ] 

Aka  js1 

X Bm  [ khm-fk-t-hm-*-!  kha-ah4-ak-a] 

Aka 

X fla  r2kh  ♦ (kfl  ) ♦ (k-1)“i 

Ak  L • J 

For  the  required  average  flow-time  of  passengers  we  must 
multiply  this  by  2.  Also,  using  (A3. 13)  we  get 


F 


k 

q;b,a 


2 

_L  fSkh  <*•  (k-»-1  ) ♦ (k-1)“| 

2lc  k^h  L ■ J 

r 

«>  q*  [2  ♦ JL(U1/k)  ♦ l(1-1/k) 

k^  L h 


[] 


This  gives  the  required  result. 
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